Towards a holographic realization of the Quarkyonic phase 
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Large- N c QCD matter at intermediate baryon density and low temperatures has been conjectured 
to be in the so-called Quarkyonic phase, i.e., to have a quark Fermi surface and on top of it a 
confined spectrum of excitations. It has been suggested that the presence of the quark Fermi 
surface leads to a homogeneous phase with restored chiral symmetry, which is unstable towards 
creating condensates breaking both the chiral and translational symmetry. Motivated by these 
exotic features, we investigate properties of cold baryonic matter in the single flavor Sakai-Sugimoto 
model searching for a holographic realization of the Quarkyonic phase. We use a simplified mean- 
field description and focus on the regime of parametrically large baryon densities, of the order of the 
square of the 't Hooft coupling, as they turn out to lead to new physical effects similar to the ones 
occurring in the Quarkyonic phase. One effect, the appearance of a particular marginally stable 
mode breaking translational invariance and linked with the presence of the Chern-Simons term in the 
flavor brane Lagrangian, is known to occur in the deconfined phase of the Sakai-Sugimoto model, but 
turns out to be absent here. The other, completely new phenomenon that we, preliminarily, study 
using strong simplifying assumptions are density-enhanced interactions of the flavor brane gauge 
field with holographically represented baryons. These seem to significantly affect the spectrum of 
vector and axial mesons and might lead to approximate chiral symmetry restoration in the lowest 
part of the spectrum, where the mesons start to qualitatively behave like collective excitations of 
the dense baryonic medium. We discuss the relevance of these effects for holographic searches of 
the Quarkyonic phase and conclude with a discussion of various subtleties involved in constructing 
a mean-field holographic description of a dense baryonic medium. 



I. INTRODUCTION 

One of the outstanding open problems in the theory of 
strong interactions is the phase structure of QCD at low 
temperatures and at intermediate baryon densities. This 
region of the phase diagram is relevant for the physics of 
neutron stars pQ and is a subject of significant research 
efforts in the heavy ion collision experiments at low and 
moderate energies [2J. 

In the vacuum and for low temperatures and densities, 
QCD is realized in the hadronic, confined phase, where 
the chiral symmetry is broken. Based on lattice QCD 
simulations [3] and 't Hooft anomaly matching [3] we 
have strong theoretical reasons to suspect that as we in- 
crease the temperature at vanishing baryon density, there 
is a single transition of a crossover nature to a deconfined, 
chirally symmetric phase - the quark-gluon plasma. 

However, little is known about the opposite region of 
zero temperature and intermediate densities, as in this 
regime of parameters it is hard to apply first princi- 
ple methods. More concretely, at moderate tempera- 
tures and baryon densities the perturbative expansion 
faces convergence issues, lattice simulations suffer from 
fermionic sign problem, and conclusions from anomaly 



* j.deboer@uva.nl 

t b.d.chowdhury@uva.nl 

•f m.p.heller@uva.nl; On leave from: National Centre for Nuclear 

Research, Hoz a 69, 00-681 Warsaw, Poland. 
§ jakubj@ift.uni.wroc.pl 



matching can no longer be drawn because Lorentz sym- 
metry is broken. 

Using qualitative large-AT c arguments, McLerran and 
Pisarski have recently conjectured that at intermediate 
baryon density and at low enough temperatures a new, 
so-called Quarkyonic, phase appears |S]. This phase has a 
Fermi sea of quarks, but the excitation spectrum consists 
of color singlets just as it does in the vacuum. Moreover, 
quarks from the Fermi sea occupy states required to form 
the chiral condensate and as a consequence, above a cer- 
tain critical density the condensate is believed to vanish 
leading to chiral symmetry restoration. This is reflected 
in the spectrum of mesons obtained using model studies 
in [51 [7], which consists of exact chiral multiplets. 

In the presence of interactions, the quark Fermi sur- 
face is unstable against the formation of quark-quark hole 
condensates [8 . This breaks the chiral symmetry again, 
but now the chiral condensate oscillates in space around 
its zero value. A preliminary analysis of this phenomenon 
in the context of Quarkyonic matter was done in [5J HOj , 
showing that it can have a significant impact on the low 
energy behavior of the system. 

A natural consequence of the existence of a new phase 
is the appearance of a Triple Point in the QCD phase 
diagram. This possibility has been explored to predict, 
both qualitatively and semi-quantitatively, the features 
of hadron production in heavy ion collisions [llj. 

During the last ten years a rapid development of gauge- 
gravity duality [T21TH] has taken place. It provides a 
novel first principle method of solving certain strongly 
coupled gauge theories in the large- N c limit. In particu- 
lar, holographic models of QCD share many features of 
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QCD in the vacuum, for example there is some similar- 
ity in the spectrum of the lowest lying mesons including 
the chiral symmetry breaking [15] . Another well-known 
feature where there is approximate agreement are the 
transport properties (like the famous ratio of the shear 
viscosity to the entropy density |16j ) in the strongly cou- 
pled quark-gluon plasma phase as observed at RHIC and 
LHC. These and other results make the gauge-gravity du- 
ality interesting not only from the point of view of QCD 
phenomenology, but also from the more fundamental per- 
spective of improving our theoretical understanding of 
non-perturbative phenomena in gauge theories (see [17] 
for a broad overview of these research trends) . 

In the gauge-gravity duality there are no fundamen- 
tal obstacles in studying holographic models at non-zero 
baryon density, or more generally in the presence of 
fermionic densities^ Hence, the results obtained using 
the gauge-gravity duality might also shed light on oth- 
erwise inaccessible quantities in cold QCD at moderate 
baryon densities, as they previously did in the case with 
the transport properties of the strongly coupled quark- 
gluon plasma. 

Our approach is based on noting that McLerran and 
Pisarski's argument about the existence of the Quarky- 
onic phase follows from large-./V c counting and so might 
be also applicable to holographic models. We hope that 
first principle investigations in string theory construc- 
tions will bring new results to the existing phenomenol- 
ogy of the Quarkyonic phase and also lead to further 
progress in understanding holographic QCD. 

There are three issues that seem to be crucial for claim- 
ing a holographic realization of the Quarkyonic phase. 
We address all of these at a preliminary level in our pa- 
per: forming a quark Fermi surface in the holographic 
model of interest when it is put at finite baryon density; 
having colorless excitations of the quark Fermi surface 
with the chiral symmetry (approximately) restored in the 
spectrum and identifying a possible mechanism leading 
to instabilities breaking translational invariance at large 
distances. 

As a holographic QCD model we choose the one in- 
troduced by Sakai and Sugimoto in [TSl HH] and study its 
properties at vanishing temperature but non-zero baryon 
density. We focus on the case of having a single quark 
flavor (Nf = 1). Our chief motivation for focusing on this 
particular realization of holographic QCD is its simplic- 
ity and the fact that it nevertheless incorporates chiral 
symmetry breaking. In the large N c limit, the Nf — 1 
analogue of the chiral symmetry breaking is the axial 
symmetry breaking and we will refer to this phenomenon 
simply as to the chiral symmetry breaking. 

The model of interest [HI [19] is based on a system of 
N c D4-branes compactified on a circle with antiperiodic 



The latter observation lies at the very heart of the AdS/CMT 
program trying to apply holography to condensed matter systems 
(see 1181 for a recent review). 



boundary conditions for fermions. Upon backreaction 
these generate a spacetime dual to the adjoint sector of 
the (3+l)-dimensional gauge theory [5U]. Nf D8-branes 
and Nf D8-branes localized at two different points of 
the circle represent Nf flavors of chiral quarks in the 
quenched approximation |21) . 

At vanishing temperature, which is what we are inter- 
ested in, the spacelike circle shrinks to zero size in the 
IR part of the geometry. The D8 and D8-branes, which 
in the UV are localized at different points on the cir- 
cle, must then connect in the IR and form a U-shaped 
configuration. This purely geometric mechanism is the 
dual description of the chiral symmetry breaking. In- 
deed, normalizable solutions of the gauge field on D8- 
branes in vacuum - corresponding to some of the spin-0 
and spin-1 mesons - have Nj massless modes interpreted 
as iV| — 1 pions and the rf meson, which in the large N c 
limit becomes the Goldstone mode of the spontaneously 
broken axial symmetry [22] . Another indication of chiral 
symmetry breaking is that the holographically obtained 
masses of the axial and vector mesons differ. 

In the Sakai-Sugimoto model, a non-zero baryon den- 
sity at vanishing temperature is realized by an antisym- 
metric configuration of the electric field in the holo- 
graphic direction on the U-shaped D8-brane. Such a 
configuration requires sources (baryons), which in the 
Abelian case considered here are D4-branes living on the 
D8-brane and wrapping the t x S 4 part of the geometry. 
After a Kaluza-Klein reduction on S 4 , the D4-branes cor- 
respond to massive, electrically charged point particles 
living in the world volume of the flavor D8-brane. 

The model studies of the Quarkyonic phase suggest 
that the chiral symmetry gets restored at sufficiently 
large baryon densities [6l [TJ (23] . On the gravity side there 
is no obvious way of achieving chiral symmetry restora- 
tion in the confining phase, as opposed to the deconfined 
phase where there is a natural mechanism to decouple 
the D8- and D8-branes by extending them all the way 
to the horizon. Therefore, if there is a dual description 
of the Quarkyonic phase it will necessarily involve some 
new ingredients. 

Although the holographic baryons are point like (sin- 
gle flavor case) or almost point-like (two flavors [23]) ob- 
jects, on the field theory side they correspond to parti- 
cles of a finite size [25] when both the number of col- 
ors N c and the 't Hooft coupling constant A are very 
large. This seems to suggest that in order to have over- 
lapping baryons on the field theory side, and hence at 
least naively form a quark Fermi surface, one only needs 
an 0(1) (in both N c and A) three-dimensional density 
of point-like holographic baryons residing at the bottom 
of the U-shaped D8-brane. Attempts to construct such 
holographic baryon densities were made in [551 12Z] • We, 
however, will focus on parametrically higher densities, of 
order of 0(A 2 ), as only in this case is there a possibility 
of having new physical effects bearing similarities to the 
ones expected to occur in the Quarkyonic phase. The 
first such effect, possible instabilities towards breaking 
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translational invariance triggered by the presence of the 
Chern-Simons term in the flavor brane Lagrangian, were 
already present in the literature on holographic QCD [2"B1 - 
130] but were never comprehensively analyzed at vanishing 
temperature and non-zero baryon density. The second, 
to the best of our knowledge, completely new effect that 
we investigate here are density-enhanced interactions of 
mesons with baryons. As we show later in the article, 
these may significantly affect the meson spectrum. 

It turns that such parametrically large baryon densi- 
ties have a peculiar effect on the distribution of the holo- 
graphic baryons (D4-branes wrapping S 4 ) on the flavor 
brane. In the discussion below we will follow |31j where it 
was noticed that due to the electromagnetic repulsion, at 
densities of order C(A 2 ), the D4-branes will fill not only 
the field theory directions, but also the holographic direc- 
tion. Such a configuration was explicitly constructed in 
[31j in a mean-field approach for the antipodal embedding 
of the D8-brane and it is the starting point for our consid- 
erations here. A further justification that the holographic 
baryon density will fill four spatial directions, including 
the holographic direction, at a large enough baryon den- 
sity on the boundary comes from the recent study [32"] . 
In this work the authors, based on large N c scaling, argue 
that the structure of the ground state of this holographic 
model of QCD is necessarily a crystal on the flavor brane. 
Such a holographic crystal structure for sufficiently large 
densities will span not only field theory, but also the holo- 
graphic directions as energetic arguments show. 

As mentioned earlier, mesons appear as normalizable 
perturbations of the flavor brane gauge field. It is natu- 
ral to expect that such modes persist at non-zero baryon 
density, albeit with possibly different properties than in 
the vacuum. We already see at this level that in the 
holographic model we can have a quark Fermi surface, 
assuming having overlapping baryons and strong inter- 
actions between their holographic counterparts is suffi- 
cient for this, and at the same time colorless excitations. 
These are the defining features of the Quarkyonic phase. 

A generic feature of holographic QCD models is the 
presence of a Chern-Simons term for the gauge field 
which represents in the bulk the chiral flavor degrees of 
freedom, and which is responsible for the flavor anomaly 
on the field theory side [T3J [33J . Recently, it has been 
shown that the presence of the Chern-Simons term to- 
gether with a non-zero value of the bulk electric field 
may generate unstable modes carrying non-zero momen- 
tum [33J [33] . The condensation of such modes leads to 
the appearance of an inhomogeneous phase |35j . Indeed, 
such an effect is present in the Sakai-Sugimoto model in 
the deconfined phase at non-zero baryon density [53] and 
also at vanishing temperature and sufficiently large axial 
chemical potential [30j . 

Although the instability pattern resembles the one 
found in the case of the Quarkyonic phase, the operators 
condensing arc different. Superficially, this seems to be 
enough to discourage linking Chern-Simons term-induced 
instabilities with the ones appearing in the Quarkyonic 



phase [HI HO]- However, it is conceivable that the instabil- 
ity due to the Chern-Simons term also triggers a modula- 
tion in the operators which are conjectured to condense 
in the Quarkyonic phase. We will not pursue this exciting 
idea here besides some brief comments in the concluding 
remarks, but we will nevertheless search for modulated 
phases due to the Chern-Simons coupling, as the break- 
ing of translational invariance will affect the properties 
of the background solution and hence the spectrum of 
mesons. We will follow [29] to search for an indication of 
an instability and for marginally stable modes, i.e., time- 
independent normalizable solutions of the equations of 
motion for the bulk gauge field carrying non-zero mo- 
mentum. It is interesting to note that the Chern-Simons 
term, which is often crucial in order to obtain instabil- 
ities, becomes relevant only for large densities, like the 
ones that we consider here, and otherwise its influence 
on the dynamics of the theory is subleading at large A. 

In this paper we take a first step in trying to under- 
stand the, so far neglected, interactions of mesons with 
the dense baryonic medium. These seem to be the sim- 
plest possible mechanism leading to a partial decoupling 
between the D8- and D8-brane in the regime of interest. 
In hindsight, the mean-field treatment of the Nf = I 
configuration adopted originally in [31] turns out to be 
a brane analogue of an electron star system introduced 
in the AdS/CMT context in [35]. The relevant action 
is that of charged dust coupled to the DBI action for 
the electromagnetic field. We argue that such an action 
captures part of the response of the baryonic medium 
to meson perturbations. The element that is not taken 
into account by the dust approximation is a restoring 
force acting on the holographic baryons once they are 
perturbed from their equilibrium position in the under- 
lying crystal lattice. 

We leave the study of this, possibly very important 
yet hard to implement in the coarse-grained description, 
effect for future investigations. We subsequently focus 
on finding zero momentum transverse fluctuations of the 
gauge field (corresponding to spin-one mesons) in the 
presence of a non-zero baryon density and show that 
the electromagnetic interactions with the D4-branes on 
the D8-brane significantly affect the frequencies (meson 
masses). In particular, the lowest mesons at large densi- 
ties turn out to be bound states in an effective potential 
generated by the presence of the baryonic charge on the 
flavor brane. This leads to the possibility that the spec- 
trum organizes itself in approximate axial (the Nf = I 
analogue of chiral) doublets. In particular, the analysis of 
this effective potential which appears in the Schrodinger 
form of the equation for in-medium meson perturbations 
suggests that such a possibility may indeed be realized 
for a particular form of the embedding, albeit only for 
the lowest axial and vector meson. We comment on the 
relevance of the restoring force that we neglected in the 
conclusions section. 

The plan of the paper is as follows. In section [TT] we 
review various features of the Sakai-Sugimoto model rel- 
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evant for our studies and explain how to incorporate 



III 



baryons in the holographic description. In section 
we show how the coupling constant (A) dependence of 
different terms in the bulk Lagrangian leads to possi- 
bly interesting physical eff ects at baryon densities that 
scale as C(A 2 ) 



In section 



IV 



we attempt to construct 
zero temperature states in the Nf — 1 Sakai-Sugimoto 
model with large baryon density for antipodal and non- 
antipodal embeddings, generalizing the results of [3T]. In 
sections [V] and |VI| we discuss interesting physical effects 
associated with large densities. In section|V]we search for 
normalizable time-independent perturbations of a gauge 
field carrying non-zero momentum, whose presence sig- 
nals an instability towards creating striped structures. 
Section |VI| is devoted to recovering part of the spectrum 
of mesons at non-zero density and explaining a possible 
mechanism of chiral symmetry restoration. The last sec- 
tion summarizes our results and discusses open problems. 



II. THE HOLOGRAPHIC MODEL 

A. Vacuum solution 

The starting point of our studies is the geometry gen- 
erated by N c D4-branes compactificd on a circle with an- 
tipcriodic boundary conditions for the (adjoint) fermions. 
Such a background, representing the ground state of the 
adjoint sector in the dual field theory, is a solution of type 
IIA supergravity with a nontrivial metric gAB, the dila- 
ton <f> and the four-form In the following we will de- 
note the non-compact D4 world volume ("field theory") 
directions by x^, the coordinate on compactification cir- 
cle of circumference St by r 



(1) 



the radial ("holographic" ) direction by u, the line element 
on the four-sphere by dfl 2 , the unit four-sphere volume 
by I/4 



(2) 



and its volume form by 64. The full solution now reads 



ds z = g AB dx A dx a = {u/R)^{n^dx^dx v + f(u)dr z ) + 

+(u/R)- 3 H^- )+ u*dnl), 



= g s (u/R) 4 and F 4 



2irN r 



V 4 



(3) 



with the curvature scale set by R, the string length and 
the string coupling denoted respectively by / s and g s and 
f(u) being 



Regularity at the tip of the cigar formed by u and r 
directions relates ukk and St 



St 



Air i? 3 / 2 

U KK 



After introducing the Kaluza-Klein mass 



2n 3 u 



Mkk = ^ = 



1/2 

KK 



St 2 i?3/ 2 : 



(5) 



(6) 



the mapping between the field theory parameters - the 
number of colors N c and 't Hooft coupling constant A — 
gy M N c - and i?, l s and g s becomes 



f 



A 3 



9* 



8tt 2 M kk R?N% 



and l 2 s = 2M KK R 3 \- X . (7) 



Note that everywhere in the text we are using the four- 
dimensional 't Hooft coupling constant A and the four- 
dimensional Yang-Mills coupling constant gym- It is 
perhaps worth clarifying here what the conditions un- 
der which we can trust the supergravity approximation 
are. As explained in detail in [37], demanding the small- 
ness of g s and l s /R while keeping MkkR fixed amounts 
to 



1 < 9ymN c < 1/9y M , 



(8) 



which in terms of the 't Hooft coupling constant A and 
the number of colors N c translates to the following set of 
inequalitie^] 



I < A 3/2 < N c 



(9) 



The key insight brought by Sakai and Sugimoto is the 
realization that, from the perspective of the dual field 
theory, adding to this background Nf coincident probe 
D8-branes and Nf coincident probe D8-branes localized 
at different points of the circular direction r at large it, 
introduces Nf chiral flavors in the quenched approxima- 
tion and also leads to the chiral symmetry breaking in the 
IR. Indeed, the form of the background geometry neces- 
sitates that asymptotically separated D8- and D8-branes 
must connect at some smaller value of radial coordinate 
and form a U-shapcd stack. In the following we will fo- 
cus on Nf = 1 case and we will consider asymptotic 
separation of D8- and D8, denoted as L < St/2, as an 
additional independent parameter specifying the model 
of holographic QCD. 

To summarize, dual field theory has two dimensionless 
parameters - N c and A (as Nf = 1) - and two scales - 
Mkk and L. Although Mkk plays the role of Aqcd 



2 In particular, in the remainder of this paper, when we discuss 
the scaling of various quantities as A — > 00 we implicitly assume 
that we keep MkkR fixed, and scale N c ~ A 3 / 2 and _R ~ A 1 / 2 
in order to stay in the supergravity regime. 
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and in particular sets the scale for glueball masses, L 
seems to have no analogue in QCD. 

The dynamics of the flavor degrees of freedom is gov- 
erned by the DBI action for, in our case, a single D8- 
brane. This object is wraps the x^ and the u directions, 
as well as the four-sphere. The degrees of freedom ap- 
pearing in the action are ten embedding functions X A 
of nine world volume coordinates a a and the U(l) gauge 
field A a which has the field strength F ab . With the in- 
duced metric given by 



9ab = gABd a X A d b X l 



(10) 



the DBI action reads 



S 



DBI 



-A*8 



J dVe" V-det (g ab + 2nV*F ah )> (11) 



where 



Ms 



(2^/9 ' 



(12) 



The action ( 11 ) needs to be supplemented by the coupling 



of the U(l) gauge field to the background C3 form given 
by 



Scs = -A* / C 3 A Tr(2< 2 F) s 

JD8 



(13) 



In the following, we will choose the gauge such that 

x a =X a (a) =a a (14) 

for all X a apart from t, which is the transverse position 
of D8-brane. The latter is given by a non-trivial scalar 
function y of world-sheet coordinates x 



(15) 



In this article we will be interested in SO(5) singlet con- 
figurations and hence in the Kaluza-Klein reduction on 
the four-sphere we will keep the zero modes only. Af- 
ter the Kaluza-Klein reduction and integrating by parts, 
Scs becomes the action for the 5-dimensional Abelian 
Chern-Simons term 



Scs = 



N c 
24ir 2 



AAFAF 1 



(16) 



where K+ and K ' denote the holographic and the field 
theory directions respectively. 

Before we move on, it is convenient to rescale world 
volume coordinates and the gauge field by 



u = Ru, x~ 
9ab = 9ab and A 



»=Rx», y = Ry, 
R 



A„ 



(17) 



Note, that rescaling A a as above and treating A a as 
an 0(1) number leads to parametrically large (in A) 
baryon and axial densities. For the latter baryon-baryon 



and baryon-meson interactions become important. Such 
rescalings were (implicitly) also made in many of the pre- 
vious works on holographic QCD at non-zero baryon den- 
sity, such as p51 l2"9l |3"T] . We will come back to this issue 
in section Mil 



Having done the rescalings (17), the five-dimensional 
action reads 



S/c=- j dud 4 xu^ 4 J- det (g ab + F ab ) 

+a [ dud 4 x e abcde A a F bc F de , (18) 



where e abc d e is the Levi-Civita symbol, the overall dimen- 
sionless factor c is given by 



c = n 8 V 4 R 9 g s \ 
or in terms of the field theory parameters by 

-AU 3 



1 



768 7 r 5 i? 3 M| rK 



(19) 



(20) 



and the coefficient of the Chern-Simons term a reads 

1 



(21) 



In the vacuum F ab vanishes and the embedding of D8- 
brane is given by r being a function of the radial coor- 
dinate only. The latter is also true at non-zero baryon 
density. As the D8-brane is shaped like the letter "U" 
and has two asymptotic regions, we will define world vol- 
ume radial coordinate z as 



u/ukk = £ ( 1 + |J 



(22) 



In that way, z running from — 00 to 00, interpolates 
between two asymptotic regions of D8-brane and £ de- 
notes the lowest radial position of D8-brane world vol- 
ume. For antipodal embedding, i.e., D8 and D8-brane 
asymptotically localized at the opposite points of the 
circle (L = St/2), symmetry dictates that £ = 1 and 
y'(z) = 0. In more general situations, £ will be a non- 
trivial function of L, determined by solving the equations 
of motion for D8-brane embedding function y{z) 

POO 

t/R = y{z) subject to L/R = 2 1 y'{z)dz. (23) 

For y being a function of z only, the induced metric on 
R 1 ' 3 x R+ reads 



g ab dx a dx b = j 



4 |7/6 2 2 



9(£2 + z 2 ) 5 / 6 (£ 3 - l+~z 2 



(24) 



^ — )\ l MM }d~ z 2 + Jae + z^^d^dx". 



G 




L/St 



FIG. 1. The plot of the lowest position of the flavor brane £ as 
a function of asymptotic separation L. The ratio L/5r = 0.5 
corresponds to the antipodal embedding. 
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FIG. 2. The relative mass difference given by (29 1 for the 
three lowest axial and vector mesons. The continuous curve 
corresponds to lowest axial and vector mesons, the dashed 
curve to the second lowest and the dotted one to third lowest. 
The difference is attributed to the axial symmetry breaking, 
which is Nf — 1 analogue of the chiral symmetry. 



The equations of motion for y{z) have the first integral 
which, for vanishing F ab , takes the form 



€ 9/4 (i + p 



s2. 19/12 



+ ^ 2 -l )y'(z) 



^-det(g ab ) 



e 3 -i (25) 



and allows us to algebraically solve for y'(z). The integra- 
tion constant on the RHS is fixed by small- z expansion 
of the LHS. Similar consideration will be used to deter- 
mine the embedding in the presence of very large baryon 
density. See Fig. [l] for a functional dependence of £ on 
the asymptotic separation L. 

In the rest of the article we will (partially) fix the gauge 
freedom in A a by demanding that Ai = leaving the \x 
components non-trivial. 

Small normalizable perturbations of the embedding 
function by and the gauge field SA^ have the interpre- 
tation of spin-zero (scalar and pseudoscalar) and spin- 
one (vector and pseudovector, here referred to as axial) 
mesons [15] . In the article we focus our attention only on 
the latter, which correspond to perturbations SA^ van- 
ishing as z — > ±oo. Such an assumption misses a sin- 
gle normalizable (in the sense of having SF a b — > as 
z — > ±oo ) mode having the asymptotics of a pure 4- 
dimensional gauge, i.e. A^ — > ±d ll (j)(x l/ ) as z — > ±oo. 
This mode is massless and corresponds to the rj meson 
being, in the large-A^ limit, the Goldstone boson of the 
spontaneously broken axial symmetry |15j . 

The masses of the vector and the axial vector mesons 
can be obtained by analyzing the spectrum of the nor- 
malizable perturbations of the form 



SA^dx* = 5h 3 {z)e- lM dx z 



(26) 



where the relations between frequencies of normalizable 
modes w and the masses of the mesons m in the dual 



field theory reads [15] 



m = -M kk uj. 



(27) 



Such an ansatz leads to mesons at zero momentum with 
polarization in one of the spacelike field theory directions, 
which without the loss of generality can be chosen to be 
Linearizing the equations of motion for the action 
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(181 leads to Schrodinger-like equation for Sa' 3 (z) and u> 



V-det g ab dz 



\ {f 9 / 12 (f + 5 2 ) 19 / 12 v/=dct^ S~a' 3 (~z)} + 



£11/12 (£2 _|_ ^2)11/1 



(28) 



where det g ab is the determinant of the induced metric 
(24 1 which encodes the dependence on the embedding. 



Imposing the normalizability condition and identifying 
the frequency of the eigenmode with the mass of corre- 
sponding meson leads to the meson spectrum in the vec- 
tor and axial vector channels. As (28) is invariant under 
z — > —z, solutions Sa 3 are either symmetric or antisym- 
metric functions of z. As explained in [15j . these corre- 
spond respectively to vector (symmetric under z — > — z) 
and axial mesons (antisymmetric under z — > —z). 

The difference in the masses of the corresponding vec- 
tor and axial mesons, i.e. lowest vector and lowest axial 
mesons, second lowest vector and second lowest axial etc, 
is attributed to the axial symmetry breaking. Fig. [2] and 
Fig. [3] show the relative mass ratios 



Am 



rriA — my 



and relative mass-squared ratios 
Am 2 = 



J m 2 , + m 2 r 



(29) 



(30) 
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FIG. 3. The relative mass squared difference given by (301 
for the three lowest axial and vector mesons. Notation is as 
in Fig. 

for the lowest three vector (with mass my) and corre- 
sponding axial mesons (with mass m^) as a function 
of asymptotic separation between D8- and D8-brane L. 
The latter quantities (squares of masses) will become im- 
portant in section |VI D[ where we analyze the equation 
for vacuum meson perturbations ( 28 ) by putting it in 



the form of the Schrodinger equation, and also consider 
its generalization to large densities. The key observation 
following from Fig. [2] is that the relative mass differ- 
ence (29) between the lowest axial and vector mesons 



is very significant, approximately 40%, and does not de- 
pend strongly on the asymptotic separation L. For higher 
order pairs of vector and axial mesons the relative mass 
difference becomes smaller and for very large quantum 
numbers becomes insignificant. This can be understood 
as the chiral symmetry restoration in the UV part of the 
spectrum [38 . 

In section |VI| we will revisit perturbations of the form 



( 26 1 in the presence of parametrically large baryon den- 



sity and investigate how some of effects generated by the 
dense baryonic medium affect the masses of vector and 
axial vector mesons. It will turn out that the effects 
considered gcnerically lower the relative mass difference 
between mesons in question ( 29 1 , in certain situations to 



as much as 10%. We will perform a thorough analysis of 



this phenomenon in section VI D 



B. Introducing baryons 

Our goal is to study the theory at non-zero baryon 
density. In the grand canonical ensemble this amounts 
to modifying its Lagrangian density by 



C ^ C - fiQ 



(31) 



where fi is the quark chemical potential and Q is the 
density of the quark charge expressed in terms of left 
and right quark operators as 



(32) 



In this formulation baryon density B can be obtained 
from the quark density Q by 



Q = N C B. 



(33) 



On the gravity side in the gauge Ai = 0, the fields dual 
to ip L ^^ipL and ipRlfj-^R are components of D8- and 
D8-brane U(l) gauge field. In other words ip L ^^ipL and 
^rI^R can be obtained respectively from z = — oo 
and z = oo asymptotics of A^iz). More concretely, the 
asymptotic solution for the time- like component of gauge 
field takes the form 



A L =C L 



and 

A |_ =C R - 



2E L 



3 V^ 

2E R 



(34) 



where Cl/ r is related to left or right quark chemical 
potential and E L / R - dimensionless electric flux on D8- 
and D8-brane respectively - is proportional to the density 
of left or right quarks q^jR in the dual field theory [3"T] 



Ql/r 



(35) 



The latter expression can be written also in terms of the 
field theory parameters as 



Ql/r 



N r X 2 



192tt a M^ k 



(36) 



As the baryon density B is a direct sum of q^ = (tp L ipL) 
and q R = (ip R ip R ) weighted by the factor 1/N C , achiev- 
ing a non-zero baryon density at a vanishing axial density 
(QA = qh — qR) requires turning on a symmetric (under 
z <H- —z) configuration of the Aq(z) gauge field. Such 
configurations give rise to an antisymmetric electric field 
which has a non-zero asymptotic flux equal to the sum 
of the fluxes on the D8- and D8-brane parts of the world 
volume. It is clear from (36) that keeping E 0(1) cor- 
responds on the field theory side to a baryon density of 
the order C(A 2 ). This is the precise meaning of the para- 
metrically large baryon densities mentioned earlier in the 
text. 

In the deconfined phase it is the bulk black hole that 
sources the flux, but at zero temperature we need to 
introduce charges (the holographic baryons) that will 
source it. Following Witten's picture of the holographic 
baryons [35]) the charges are N c end points of strings 
stretched between the flavor D8-brane and a D4-brane 
wrapping 5 4 part of the geometry. The number N c fol- 
lows from the Gauss law on S 4 , including the flux gener- 
ated by the stack of N c D4-branes, which give rise to the 
adjoint sector. 



The tension of the D4-D8 strings will cause the bary- 
onic D4-branes to reside within the world volume of the 
D8-brane 40J . Hence, from the point of view of the 
holographic and field theory directions, the baryonic D4- 
branes appear as point-like particles carrying N c units of 
charge minimally coupled to the U(l) gauge field 



Nr. 



A a —dT. 
CLT 



(37) 



In order to obtain the mass, consider the baryon world 
line with t being the parameter on the geodesic. The 
action for this is identified with the DBI action for a D4- 
brane. There is a subtle issue here that the D4-branes are 
wrapping S 4 whose size depends on the radial position. 
This leads to a mass which is a non-trivial function of 
the radial coordinate. The DBI action for the D4-brane 



Sd4 — — 



V4R 4 «4 



dW-gtt, 



(38) 



where we are keeping yj—gtt under the integral in order to 
compare with the general action for a point-like particle. 
After taking into account the fact that the action for a 
point-like massive particle at rest (i.e. it* 
the world line parametrized by t is 



J=) with 



dt (mV-fftt) , 



(39) 



the mass and the charge of the holographic baryon read 



R-fu{t)i 
6ir 2 g s l 5 s 



= N e , 



(40) 



Note, that in order to interpret m from (40) with u(t) 



ukk as the mass of the baryon on the field theory side for 
the antipodal embedding of the D8-branes, one needs to 
take into the account the redshift factor \J—gtt evaluated 
at u = ukk- This leads to an expression for the baryon 
mass matching the ones in [3T] and [23] (the latter in the 
leading order in the large A expansion). The discussion 
generalizes trivially to non-antipodal embeddings. 



III. MOTIVATIONS FOR GOING TO 
PARAMETRIC ALLY LARGE DENSITIES 

In this section we provide the rationale for studying 
holographic QCD at parametrically large baryon densi- 
ties i.e. of order of 0(A 2 ). It turns out that there are 
several reasons for doing so. 

In the first place, as explained in the previous section, 
0(A 2 ) baryon densities are dual to O(A) gauge field on 
the flavor brane. Such parametrically large gauge fields 
are the biggest ones that can be achieved on the flavor 
brane in the DBI approximation and as such require tak- 
ing into account the nonlinear nature of the DBI action. 



This means, in particular, that the equations of motion 
for mesons (and also for any other perturbations) will 
be modified leading to possibly different in-medium me- 
son spectrum. Such an effect was preliminarily studied 
in |27j with an explicit simplifying assumption about the 
charge distribution on the flavor brane, namely that all 
the charge will be localized at z = 0. We will reconsider 
this effect in section |VI[ taking into account dynamical 
charge distributions. 



It is also easy to notice, by examining the action (11) 
(13), that only for 0(A 2 ) baryon densities the Chern- 



Simons term contributes to the equations of motion on an 
equal footing with the kinetic term. This has profound 
consequences for the low energy physics of cold and dense 
holographic QCD, as the Chern-Simons term in the pres- 
ence of an electric field is known to lead to instabilities 
towards creating striped phases [291 ETH [Ml [35]. Such 
instabilities modify the background and lead to differ- 
ent properties of the excitations. Although the operators 
condensing are different from those that condense in the 
Quarkyonic chiral spirals [S], it is conceivable that the 
latter also condense in the expected way, but this effect 
might not be computable using the observables captured 
by the DBI + Chern-Simons action alone. Some further 
discussion of this point is presented in the conclusions. 
We will investigate the existence of instabilities by ex- 
amining the presence of the marginally stable modes in 
section [V] 

On a related note, the Chern-Simons term leads to 
a particular mixing pattern between vector and axial 
mesons at non-zero momentum in the presence of back- 
ground baryon density [33] . Again, this can only happen 
for 0(A 2 ) densities. We will return to this phenomenon 
in the conclusions. 

The most striking feature of our model at 0(A 2 ) 
baryon densities however is the way the holographic 
baryons organize themselves on the flavor brane. In- 
stead of forming three-dimensional structure residing at 
z = 0, very strong repulsion between the holographic 
baryons leads then to four-dimensional volume on the fla- 
vor brane, which is occupied by the holographic baryons. 
This was noticed already in [31], where the authors con- 
structed holographic background at 0(A 2 ) baryon densi- 
ties postulating a particular mean field description. We 
will reintroduce this model in the next section and in 
the remaining part of this section we will focus on impli- 
cations of having parametrically large baryon densities, 
which follow from more microscopic considerations. 

The elements of microscopic understanding came 
about very recently in |32j , where the authors considered 
simplistic models of baryonic lattices and preliminarily 
showed that for the same baryon density a lattice extend- 
ing into holographic direction has lower free energy than 
lattices residing at the bottom of the flavor brane. The 
result of these studies are charge distributions formed out 
of a very few layers residing on top of each other in the 
holographic direction. The work of |31) needs to be un- 
derstood as postulated continuum model of the situation 
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in which very many such layers are formed. 

The crucial observation behind our work is that for 
0(A 2 ) baryon densities the interactions of mesons with 
baryonic medium are enhanced to the level that they 
influence the spectrum of mesons. Although a single 
baryon's response to the flavor brane gauge field pertur- 
bation is A-suppressed due to the factor of A sitting in its 
mass, large number of baryons residing in a fixed holo- 
graphic volume makes the overall contribution an 0(1) 
effect, a priori equally important as the kinetic term for 
the gauge field perturbations. We will describe this phe- 
nomenon in detail in section IVII 

We want to stress at this point that currently there 
exists no fully reliable macroscopic description of 0(X 2 ) 
baryon densities. Because of this, we will reuse the 
coarse-grained description of 31 . In the next section we 
review the idea behind it and point out possible caveats. 

IV. ATTEMPT TOWARDS CONSTRUCTING 
MEAN FIELD DESCRIPTION OF DENSE HQCD 

A. Mean-field Lagrangian 

Having introduced the holographic baryons into the 
model and understood the motivation for studying para- 
metrically large 0(A 2 ) baryon densities, we can proceed 
to a corse-grained description of the bulk physics. 

A single baryon sits at the bottom of the D8-brane, 
as this minimizes its gravitational potential energy com- 
ing from the D8-brane being non-trivially curved. As 
more baryons are added, initially they all stay at the 
bottom of the D8-brane and get distributed only along 
the field theory direction^ as climbing the holographic 
direction costs A-enhanced gravitational potential energy. 
However, for sufficiently large densities, strong repulsion 
between the holographic baryons will start organizing 
them not in three-, but rather four-dimensional struc- 
tures [3T]l3"2]. 

This can be understood in simple terms by noting that 
the baryons residing at the bottom of the D8-brane gen- 
erate electric field with a non-zero radial component ori- 
ented towards the boundary. Imagine now adding a single 
baryon on top of the others already residing at z = and 
distributed in the field theory directions. If the density 
is large enough, then the electric repulsion between the 
z = baryons and a new baryon that we are trying to 
add to the system is such that the equilibrium position 
for the added baryon is no longer at z = 0, but rather at 
some \z\ > 0. As the baryon mass scales as N C X and its 
charge as N c , representing non-zero baryon density in the 
dual field theory in terms of the three dimensional den- 
sity of the holographic baryons residing at z = 0, which is 



3 As in the Ne = 1 Sakai-Sugimoto model baryons repel each 
other, one should think about the dual field theory at finite 
baryon density as being defined in finite volume, e.g. on a torus. 



the description adopted in the references [26, 27J, will be 
the relevant way of describing the system only for 0(1) 
baryon density in the dual field theory. 

For parametrically larger values of the baryon density, 
the bulk description in terms of four-dimensional density 
of the holographic baryons will be appropriate. Such a 
way of describing the system was adopted for the first 
time for antipodal D8-brane embedding in the reference 
[31] . and here we will generalize this picture to the non- 
antipodal case. 

In the reference [3T] it was also conjectured that the 
holographic baryon density being not three- but rather 
four-dimensional is a dual manifestation of building up 
the quark Fermi surface or even the direct gravity dual 
of the Fermi surface. Indeed, the bulk charge distribu- 
tion climbing up in the bulk radial direction, which is 
associated with the energy scale in the dual field theory, 
makes the connection with building up the Fermi surface 
tempting. There is, however, another motivation for this 
association. Although the baryon size (on the bound- 
ary) is 0(1) in both N c and A, and hence baryons can 
significantly overlap at relatively low, three-dimensional 
bulk densities, for such configurations the holographic 
system is not in a strongly correlated phase. The lat- 
ter feature, strong interactions between the holographic 
baryons, arises precisely at 0(A 2 ) densities, which is in 
line with the point of view presented in |32j . 

The holographic system corresponding to parametri- 
cally large baryon densities, after the Kaluza-Klein re- 
duction on the four-sphere, is a set of large number 
of charged point-like particles interacting via the U(l) 
gauge field with DBI kinetic term. We are interested in 
the coarse-grained description of this system, as superfi- 
cially only long- wavelength features will be important for 
understanding instabilities associated with the presence 
of the Chern-Simons term and non-zero density, as well 
as meson perturbation whose frequencies in the vacuum 
are 0(1) both in N c and A. 

There are general arguments following from large N c 
counting, that the ground state structure of such a sys- 
tem is crystalline [35] . These arguments are based on the 
kinetic energy being suppressed by the iV c -heavy mass of 
a single (holographic) baryon while the potential energy 
is enhanced by the 7V c -large quark charge. Below we will 
review the analysis of [3T], where, setting aside certain 
subtleties that we will discuss, the macroscopic baryon 
densities - four-dimensional holographic crystals - were 
constructed in a simple coarse-grained description. The 
key idea there was to minimize the total energy of the 
dynamical electromagnetic field generated by the contin- 
uous charge distribution, as well as the potential energy 
of the charges in the dynamical electromagnetic and ex- 
ternal gravitational fields. Such a description bears sim- 
ilarity to the AdS/CMT electron star construction [36] 
and in some respects can be understood as a probe brane 
analogue of it. 

In the coarse-grained description, the relevant degrees 
of freedom will be the gauge field (as it encodes the axial 
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and the vector currents), as well as the quantities de- 
scribing a continuous charge distribution, i.e. the local 
energy density e and local charge density p. Note that 
the gauge field we will be talking about will be a mean- 
field one. Nevertheless we will explicitly assume that its 
kinetic term is still given by the DBI Lagrangian, which 
is the case in the vacuum. 

As the system is put in an external gravitational field 
(the curved D8-brane world volume), we will take into 
account the redshift by introducing the local velocity of 
the charge density u a , which is normalized as 



and (ii) assuming that we are in flat space. An approxi- 



1. 



(41) 



Although u° encodes the red shift factor which is of or- 
der 0(1) in N c and A, it is important to note that we 
will never be interested in having 0(1) components of 
the velocity field in the radial direction, as well as in the 
spacelike field theory directions (up to boosts, which are 
trivial). The reason for this is that the crystal- like na- 
ture of the configurations we study would not allow for 
such motion. However, infinitesimal, of order 0(1/A), 
displacements in spacelike directions are not only per- 
mitted, but also lead to new interesting physical effects 
that we will preliminarily study in section fvT} 

The missing ingredient in our description is the equa- 
tion of state that relates the energy and charge density. 
For a free system both quantities are proportional to the 
number density with the proportionality constants being 
respectively the mass and the charge of a single parti- 
cle. In the interacting system, the energy density can be 
corrected by interactions between the constituents in a 
given infinitesimal volume. Another way of looking at 
this is that by using the coarse-grained gauge field in the 
mean-field approach we are no longer including in the 
total energy count of the system highly oscillatory be- 
havior of the microscopic gauge field coming from gran- 
ular structure of the substance under studies. In order 
to investigate this issue better, let us assume that the 
local energy density e and the local charge density p are 
proportional to the local baryon density n 



mb(z)n and p = q b n. 



(42) 



and estimate the mistake made by neglecting interactions 
within the small volume element SV . 

Note that only for a holographic number density of 
order 0(A 2 ), the energy density, the potential energy of 
charges proportional to pA and the kinetic term for the 
gauge field all contribute on equal footing. The 0(A 2 ) 
density of the holographic baryons results in a holo- 
graphic lattice spacing d of order 0(A -1 / 2 ). We now 
wish to examine whether the mean field description is a 
reasonable approximation. To do so, we should compare 
the interactions of a collection of lattice points in a small 
volume SV to those of a smeared charge distribution. 

The interactions of a collection of lattice points can be 
estimated by (i) assuming the interactions are not those 
of a DBI theory but of five-dimensional electromagnetism 



V.**; """"""6 """" 

mate expression is thus 



i n 



sv 



A d 2 ( P 2 



-Pi 



-Pi 



-Pir 



(43) 



where the pj 's denote integer node numbers in the z and 
x 1 directions and we sum over all the holographic baryons 
residing in the volume SV. The sum over all nodes in (43 ) 



diverges as the square root of the total number of nodes 

Anodes 



Se « N c Tly/ # no des, 

or after using Anodes ~ nSV as 

Se « J\U 3/ W 2 . 



(44) 



(45) 



If we replace the lattice by a continuum distribution, we 
will get an answer that scales in a similar way, and so 
will the difference between the two. There are two ways 
of looking at (44 1 and (45). As n is 0(A 2 ), for any finite 
SV the changes in the energy density given by ( 45 1 are 



comparable to the energy density we started with, i.e. 
0(N C \ 3 ). Hence, if we are supposed to evaluate (451 for 
some finite/fixed volume SV, then the equation of state 
will receive order one corrections and the mean field ap- 
proximation should not be used. On the other hand, and 
this is the approach we take in the remaining part of the 
text, if we are supposed to set the volume to be infinitcs- 
imally small, then the difference in energies ( 45 ) becomes 



negligible. An alternate way of phrasing this is that for a 
large but fixed number of nodes the energy in ( 44 ) will al- 



ways be coupling constant suppressed with respect to the 
energy density of the free system. Note finally, that the 
volume-divergent contribution of the kind ( 43 ) would be 



present for arbitrarily weak electromagnetic interactions, 
however the weakness of the interactions then guaran- 
tees that evaluating it for a SV of order 0(1) still leads 
to a negligible contribution. This is precisely the rational 
adopted in the electron star construction 36J. 

Having discussed a subtlety related to the equation of 
state, we will stick with the equation of state ( [42] ) and 
in the remaining part of this subsection we will discuss 
the action principle describing our system. As the holo- 
graphic baryon density n is by definition non-negative, it 
will be useful, in the following, to express it as 

n = w 2 . (46) 

Keeping in mind that we will never be interested in ve- 
locities u a having finite components in the radial and 
spacelike field theory direction, the action principle that 
describes our system is that of a charged dust [41] 



Sdust = J d b xy / -&et(g ab ){ ~ m b {z)w 2 + 

+q b w 2 u a (A a -d a 4>) + X(u a u a + 1) } , (47) 

where q b w 2 u a A a is the standard coupling of the charge 
current q b w 2 u a to the U(l) gauge field A a , <p transforms 
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as a phase under U(l) gauge rotations and ensures on- 
shell conservation of the current, whereas the introduc- 
tion of A takes care of the on-shell velocity normalization 
(41). After introducing rescaled variables (17), the dust 



action becomes 



1/12 _ 2 
W 



S d ust/c = j d 5 xy / -det(g ab ){ - ^| 1/4 (1 + y) 11 

+ 7 w 2 u a (A a -d a ^j + \(u a u a + 1) } , (48) 



where 



2nl 2 

13 = An 2 . 7 = 12tt 2 and w= — £w. (49) 

R z 

The rescaling of other variables will not be important for 
our results, but can be easily worked out. Note again 
that the rescalings (49) imply that the 0(1) rescaled 
baryon number density w 2 corresponds to an order 0(A 2 ) 
physical baryon density w 2 . As anticipated earlier in the 
text, rescalings (fi~7| and (49) take us directly to the very 



dense phase with baryons climbing up the holographic 
direction. For rescaled densities being A-suppressed, the 
terms describing baryon dynamics drop out and baryons 
then reside at z = and appear as boundary conditions 
for the flux, exactly as in the references [551 H3- We 
think however that in the latter cases the backreaction 
of the baryons on the embedding should be treated as a 
A-suppressed effect. Otherwise, instead of taking it into 
account, one should rather pass to the effective descrip- 
tion of four-dimensional holographic baryon densities, as 
being done here or originally in the reference [31] . 

The full action is the sum of Sdbi arL d Sdust an d its 
variation leads to the equations of motion 



80 ' 

Sw 
5_S_ 

~SX ' 

SS 

5u a 



Vahw 2 U a ) = 0, 



(50a) 



/?e 1/4 (l + = !wu a {A a - dj), (50b) 



r yw 2 



(4, 



2Xu n = 0. 



(50c) 
(50d) 



These equations need to be supplemented with the DBI 
generalization of Maxwell's equations for A a 



d 



d 



1/4 



:i+£}"' 2 



-dct(.g, 



dx- d(d a A b 
= -V-detfeh^V -3ae acd efF cd F e f 



cd 



F, 



cd , 



(51) 



which comes from varying the full action with respect to 
A a , and the equation of motion for y entering the action 
through the induced metric g ab . 



Note that equations (50d) and (51) imply that there 



is an explicit coupling between the gauge field, in partic- 
ular the gauge field perturbations corresponding to the 
mesons, and the holographic baryon density. The sub- 
tlety that the velocities in the spacelike directions cannot 



be taken finite due to the underlying crystal structure is 



explained in section VI There might also be gradients 
that we neglect at this point and we elaborate on their 
role for meson perturbations in section |VI| 

Setting aside those subtleties, in the next subsection 
we will use the action (48) to construct homogeneous 



large baryon density states for the family of models 
parametrized by the asymptotic separation L of D8- and 
D8-brane. We will show that the resulting background 
(in the antipodal case) is exactly the one studied in |31) . 



B. Solutions for embeddings 

In equilibrium, the most natural assumption for the 
long- wavelength description of a non-zero baryon density 
state is that of isotropy and homogeneitjj^J This means 
that none of the functions appearing in the action ( 18 ) 



(48) depends on field theory directions x M . Moreover, 



in the rest frame, the symmetries dictate that the only 
non-trivial component of the gauge field is Aq and the 
same is true for the baryon velocity, i.e. u . 

The velocity needs to be subject to the normalization 
condition ( 41 ) , which in this situation fixes its form com- 
pletely 



£l/4(£2 + 22)1/4 



(52) 



Substituting this result into the action (47) and setting 
(f> to results in an effective action for the embedding 
function y, the baryon number density w 2 and the gauge 
field Aq of the form 



S, 



;///c =|^{-f/ 4 (i+£) 1/1 



e 

^-dci(g ab )-e' 2 {l + tf /2 A' {~ Z y- 

PV-dot(g ab )i^(l + t) 1/12 w(S) 2 + 
S 

7V-det(.g Qb )r 1/4 (l + £)- 1/4 ^) 2 i (z)}, (53) 

where det(g a f,) is the determinant of the induced metric 
(24). The action (53), in the antipodal case (i.e. y'(z) = 
Ojjis identical to the one used in the reference [3T] after 
suitable redefinitions of what is meant by the mass of the 
holographic baryons and their densitjj^] 



4 We will try to relax those in the next section, when looking for 
marginally stable inhomogeneous modes. 

5 Compared to I31| . in our approach the mass of baryonic D4- 
brane comes without the factor of y 7- 900 and u° = 1/ y/— goo is 
not absorbed in what is meant by the density (u> 2 ). Moreover, 
in our effective Lagrangian the electromagnetic coupling comes 
with u° = goo factored out. 
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We are interested in studying the theory at non-zero 
baryon density and vanishing axial density. We thus fo- 
cus on the symmetric (with respect to z f> -z) configu- 
rations of Aq(z), which, through flux conservation, need 
to be supported by a (symmetric) profile of w{z). 

The curvature of the D8-brane will make the holo- 
graphic baryons fall towards z — and the U(l) inter- 
actions between them will counter this effect. So the 
expectation is that the holographic baryon density will 
reside at the bottom of the U-shaped D8-brane. In the 
following, we will derive the properties of the holographic 
charge distribution from the equations of motion for the 
effective action (18). A technical improvement, com- 
pared to |3lj^) is the use of positivity of number den- 
sity (n(z) = w(z) 2 ) to explicitly extremize the action. 
We will do this by solving the equations of motion sepa- 
rately inside and outside holographic baryon density and 
matching the solutions. 

Indeed, let us start by considering the equation of mo- 
tion obtained by taking the functional derivative of ( 18 ) 



with respect to w(z), which after trivial simplifications 
reads 



= 0. 



(54) 



If w(z) = 0, this equation is trivially satisfied. However, 
for w > (i.e. in the presence of macroscopic density of 



the holographic baryons) (54) fixes the form of Aq{z) 

(55) 



Z 2 , 1/3 



The ratio of j3 to 7 is 1/3, but in (55 ) we will keep /3 and 
7 as explicit parameters. Our motivation for doing this 
is that it is conceivable that a more precise microscopic 
description of the system will lead to an effective mean- 
field description not dissimilar from the one we have been 
employing, but with different values for /? and 7, and as 
we will see this may lead to important qualitative changes 
in the behavior of the system. 

Consider now the equation of motion for y{z). This 
quantity appears in the effective action ( 18 1 only through 
the determinant of the metric. Note also that the matter 



part of (18) vanishes on-shell by (54). This implies that 



the embedding function y{z) knows about the matter on 
the D8-brane only through the electric field sourced by 
it. 

The equation of motion for y{z) again has a first inte- 
gral, which reads 

i 9/4 (i+f) 19/12 (f 3 +^ 2 -i) m 



-det(g ab )- £9/2(1 + ff /2 A' (~zY 



I 3 - I- (56) 



Apart from considering also non-antipodal D8-brane embed- 
dings. 



The RHS is determined from the power series expansion 
of the LHS around z = using ( 55 ) . The equation ( 56 ) 



allows us to algebraically solve for y'(z) in terms of £ and 
A' (z). 

The only remaining equation of motion is the one for 
A (z). After using (56), it reads 



d 



~2 dS dA' a {z 




\ £10/3 _ £19/3 
7w(5) 2 



5/3 



4z 2 (i 2 + - 9£ 5 / 3 (£ 3 - l) A' Q (S) 2 



5 /3 

^ £10/3 _ £19/3 + (£2 + ^\ fja + |~a _ {j 



(57) 



Inside the volume occupied by holographic charges this 
equation specifies the charge density as the gauge field is 
given by (55), whereas on the outside, the RHS is zero 



and the equation has a first integral being the dimension- 
less flux E through one end of the D8-brane 



d 



dA' (l 



1~, 
3 



|5/6 ^2 + 



3/2 



4 £5/3 5 2 _ 9^/|2 +5 2 (£3 + _ i^j A' {z) 2 



(58) 



= E, 



\ 1 10 / 3 - 1 19 / 3 + (e 2 + 5 2 ) (e 3 + is 2 - 1) 

or equivalently the baryon density B in the dual field 



theory (32) 



D = 



9&-K i M 2 (K R^ 



E. 



(59) 



The full solution is now obtained by matching the elec- 
tric field, i.e. A' (z), inside the charge distribution on the 
D8-brane with the electric field of constant flux outside 
it. The UV parameters are L — the asymptotic separation 
of the D8 and D8-brane - specifying the dual theory and 
E specifying the state. The chemical potential in this ap- 
proach is a secondary object, as it arises from the integral 
of A'q(z) from z = to z = 00 with the boundary con- 



dition that A o (0) is given by the formula (55) evaluated 
at z = 0. There are two IR parameters in the procedure: 
Zq - the radial position on the D8-brane where charge 
density terminates - and £ - the radial position in the 
target space of the lowest point of the D8-brane. 

The solutions parametrized by L and p (or equivalently 
the chemical potential fl) are obtained numerically by 
specifying p and searching for £, such that the formula 
(23) with y'{z) given by (56) is satisfied. Note that for 
each checked £ one still needs to match the electric field 
at the boundary of the holographic charge distribution. 
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FIG. 4. The electric field A' (z) as a function of the radial 
coordinate z in the antipodal case (£ = 1) for three differ- 
ent rescaled baryon densities: E = 1 (dotdashed), E — 100 
(dotted) and E — 500 (dashed). The envelope (solid curve) is 
the electric field inside the charge distribution. The disconti- 
nuity in the derivative corresponds to the radial position at 
which the holographic charge density terminates. The electric 
field never gets large, as matching condition bounds it from 
above. Note that the volume occupied by holographic charges 
increases with the baryon density on the field theory side, but 
the distribution of the charge in the core is not altered by the 
outer layers, which is why we used the envelope. 



Fig. [4] shows the profile of the electric field for antipo- 
dal embedding for four different charge densities, whereas 
Fig. [Hjdepicts the same quantity for E = 500 and four dif- 
ferent embeddings from antipodal to very non-antipodal. 
What one can easily notice is that the electric field is 
bounded from above and never gets large enough to suc- 
cessfully compete with the induced metric in the DBI 
action. This can be easily understood, as the electric 
field is continuously matched from the electric fields in- 
side and outside the charge distribution. In such a case, 
the electric field is bounded by the (smaller of) the max- 
imum value(s) of the electric field inside or outside the 
charge distribution. By looking at Fig. [4] and [5] we see 
that the relevant number is about 0.15 corresponding to 
the point where the electric field inside the charge dis- 
tribution has a maximum with vanishing first derivative. 
Indeed, although the position at which the second deriva- 
tive of ( 55 ) vanishes depends on the embedding through 
£, the value of the electric field at the maximum is fixed 
and reads 



2 1 /3 v /3 7 



0.153. 



(60) 



As argued earlier, the embedding feels the charges on 
D8-brane only through the DBI term, hence it is never 
gets significantly distorted from its vacuum form by the 
presence of baryon charges. The smallness of the electric 
field also has profound consequences on the stability of 
the system, as discussed in the next section. 

Although for the top-down model of Sakai and Sugi- 



FIG. 5. The electric field A' (z) as a function of the radial co- 
ordinate z at fixed baryon density [E = 500) for four different 
embeddings ranging from the antipodal (L/St — 0.5) to very 
non-antipodal (L/St = 0.05): L/St = 0.5 (solid), L/St = 0.2 
(dotted), L/St = 0.1 (dotdashed) and L/St = 0.05 (dashed). 
As in Fig. |4j a discontinuity in the derivative corresponds to 
the radial position at which the holographic charge density 
terminates. Also in this case the electric field never gets large 
due to the matching condition at the boundary of the bulk 
charge distribution. 



moto the values of j3 and 7 are fixed and are given by 
(49), one can also treat them in the bottom-up spirit 
as free parameters. In particular, such an approach al- 



lows us to obtain larger values of the electric field on the 
D8-brane than 0.153, and at the same time denser distri- 
butions of the holographic baryons. This has important 
consequences for the stability of the system, as well as the 
form of the meson spectrum at non-zero density. Note 
also that for completely arbitrary /3 and 7 we are not 
guaranteed that the matching condition (i.e. demanding 
continuity of the electric field) at the bulk boundary of 
the holographic charge distribution will be satisfied. 

Before we conclude this section, let us comment on and 
stress some aspects of the physics described by the model 
considered. 

In the first place, note that the gauge field inside the 
volume of charges depends mildly on the total amount of 
charge - the latter appears through the matching condi- 
tion at the boundary of the charge distribution z = zq 
constrained by the fixed asymptotic separation of the D8- 
brane. This, through equation (57), implies that the way 
the charge is distributed does not change significantlj|j 
with the increase of the baryon density B. Instead the 
charge added to the system accumulates on top of the 
existing one without changing much the original distri- 
bution. Although naive intuition would suggest that the 
charge at the bottom gets more and more compressed as 
we add more of it to the system, it is not correct and does 



7 For antipodal embedding it does not change at all, for non- 
antipodal ones it changes through the change of £ for each B. 
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not take into account the fact that there are also charges 
on z < part of the U-shaped D8-brane countering this 
effect. 

Let us consider a simple toy model which illustrates 
this effect: a system of charged point particles in four 
spatial dimensions with 1/r 2 interactions, and with a 
harmonic oscillator potential in one of the four space- 
like coordinates which we call u. One can easily show 
that in this toy model all finite density configurations 
have exactly the same constant density and are only dis- 
tinguished by their extent in the M-direction, similar to 
what we see in the Sakai-Sugimoto model. 

It is also very interesting and puzzling that in the bulk 
we are describing a finite density system with fermionic 
constituents (quarks) in terms of geometric objects (D4- 
branes) that in the vacuum correspond to a bound states 
of the fermions. Note that from the point of view of 
the dual field theory baryons are (heavily) overlapping, 
nevertheless in the bulk this seems to be taken into ac- 
count by a reasonably simple course grained description. 
The caveat might be that such a description is simple 
because it misses some ingredients. For example, in the 
non-Abelian case with Nf > 1 the holographic baryons 
are described by instantons on the stack of the flavor 
branes and these presumably interact through compli- 
cated many body interactions stemming from the non- 
linear nature of the DBI action. 



V. HOMOGENEITY OF THE DENSE HQCD AT 
LARGE DISTANCES 

A. The role of the Chern-Simons term 

So far in our description of cold baryonic matter in the 
Sakai-Sugimoto model we explicitly assumed homogene- 
ity at large distances both in the charge distribution, as 
well as in the ansatz for the mean-field gauge field. These 
usually are sensible assumptions to be made, as typically 
gradient terms increase the (free) energy. However, it was 
recently found that the presence of the Chern-Simons 
term in holographic Lagrangians for U{1) gauge fields, 
like in our case, together with large-enough radial elec- 
tric field may lead to the appearance of tachyonic modes 
of the gauge field carrying non-zero momentum |34j (see 
also [33J for earlier observation of this sort). The conden- 
sation of such modes, after including nonlinear effects, is 
supposed to give a new inhomogeneous and stable (or at 
least having lower free energy than the homogeneous one) 
state. This was indeed verified on particular examples in 
[35] and [2T5]- 

The effect in question seems to be generic and was 
found within a class of bottom-up holographic models 
having as a solution the Reissner-Nordstrom black hole 
in AdSs [34] and was also shown to exist in the very same 
setting that we are considering - the Sakai-Sugimoto 
model, albeit in the deconfined phase at large enough 
baryon density |29j . as well as at zero temperature and 



big enough axial chemical potential [3Qj. It also turned 
out that similar phenomenon can be observed in four- 
dimensional holographic models including pseudoscalar 
fields and parity- violating terms for the gauge field |42j . 
Subsequent developments in this active area of research 
include [43H48]. 

We are not going to perform comprehensive analysis 
of linear stability of our effective model, as this has not 
been achieved yet in holographic systems of less complex- 
ity than ours. Instead, following earlier works on the role 
of the Chern-Simons term in the Sakai-Sugimoto model 
[2U1 130] . we will search for the presence of a particular 
normalizable mode, such that it carries non-zero momen- 
tum, has a zero frequency and its polarization ensures de- 
coupling from other perturbations at the linearized level. 
Such a mode is called a marginally stable one, as typ- 
ically perturbations carrying momentum evolve or dis- 
perse, whereas this one stays constant in time. It is ex- 
pected^] that in the direct neighborhood of parameters 
for which marginally stable modes appear there will be 
tachyonic modes. In that sense, searching for marginally 
stable modes is testing some aspects of the stability of 
the system towards creating structures breaking transla- 
tional invariance. 

Based on previous studies in [221 |3Q], we will search for 
the presence of the following marginally stable mode 



SA^dx^ = h(z) ^cos(kxs)dx 



sin(fci 3 )di 2 | , (61) 



which will be a small perturbation on top of the finite 
density background derived in section [IV B| Such a nor- 
malizable mode turned out to be present in the Sakai- 
Sugimoto model in the deconfined phase for large enough 
baryon density and in the confined phase for large enough 
axial density. In both situations it was necessary to con- 
sider parametrically large, of order C(A 2 ), densities, as 



expected on general grounds and discussed in section III 



It is worth stressing already at this point that it was 
not the large rescaled density of charge which was really 
necessary in both cases, but rather the large bulk elec- 
tric field associated with it. Note however, that in our 
case the electric field is bounded from above by a rather 
small value (60). The second comment in place is that 



the ansatz invoked (61 ) leads to the appearance of a non- 



zero magnetic field, but it does not modify the electric 
field in any way. As the magnetic field couples only to 
the moving charge, the charge distribution is not going 
to be altered by the mode (61 ). 



The breakdown of translational invariance at non-zero 
baryon density due to the presence of the Chern-Simons 
term in the flavor brane Lagrangian might not be an 
unwanted feature after all, as the Quarkyonic phase is 



8 We verified this in the Sakai-Sugimoto model put at non-zero 
axial chemical potential and expect it to be correct in other cases 
too. 
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expected to be inhomogeneous too. Although the con- 
densate there [9] which was given by 

(■4;exp(2inx 3 joj 3 )tp) (62) 

is clearly different from the prospective one in our case 



(ip 7i cos(fci 3 ) + 72 sin(fci 3 ) ip) 



(63) 



it is not entirely unconceivable that the modulation in the 
vector or axial current triggers the one encountered in the 
Quarkyonic phase. We will come back to this issue in the 
conclusions. Last but not least, it is also worth stressing 
that any inhomogeneities at large distances will influence 
the spectrum of mesons and as the spectrum is another 
feature we investigate, the presence of instabilities is of 
clear importance and interest for us. 

To clarify the mechanism for the instability, we focus 
on the antipodal case where the expressions are less clut- 
tered. Searching for the marginally stable mode (61) in 



the non-antipodal case is a simple generalization of the 
antipodal situation, hence we skip the technical details 
and just present the results. It is instructive to derive 



the reduced action for perturbations by (61 1 by expand 



ing our coarse-grained action to the quadratic order in 
h. As expected on physical grounds, we find that pertur- 
bations in question decouple leading to the reduced Id 
Lagrangian density 

C = a{~z)~ti{z) 2 + b(S)h(z) 2 + 12a k ~h{zf A' Q (z) , (64) 



where 



a(z) = -- 



3(1 + z 2 ) 



b(z) 



2JA-9(l + z 2 y/ 3 A' (z) 2 



9(l + z 2 ) 1 /3i' Q (5)2 



6(1 + Z 2 )3 



(65) 



(66) 



and Aq(z) is the background gauge field. It is easy to 
see already at the level of equation ( 64 ) that the Chern- 
Simons term - the one entering with a coefficient - can 
lead to tachyonic behavior at non-zero momenta: al- 
though b(z) playing the role of the mass term is negative 
or zero as it needs to be, the Chern-Simons term contri- 
bution does not come with a definite sign and in certain 
circumstances can win over the standard mass term. 



Up to this point, while deriving (64 1, we did not as- 



sume anything about the form of the background gauge 
field - it could come from the axial or baryon charge or 
both. Let us now discuss a subtlety that has not been 
present in the axial case analyzed in [30] and complicates 
studies of the baryonic case here. It is easy to spot that 
a (z) and b(z) are even under parity, i.e. z — > — z. We 
see that if A Q (z) is odd under parity, which corresponds 
to purely axial chemical potential, the Lagrangian (64) is 



even under parity. This implies that solutions h(z) have 
definite parity, hence it is sufficient to construct them for 



z > 0. However, for any other form of the background ra- 
dial electric field, in particular for the vector form of the 
potential we are interested in, the solutions h{z) will not 
have a definite parity. This implies that while searching 
for normalizable solutions it is necessary to check normal- 
izability at both boundaries: z = — oo and z — oo and 
complicates searching for marginally stable modes in our 
case. The discussed feature was not taken into account 
in earlier studies in [3?5] , where the authors attempted to 
check the stability of finite baryon density configuration 
of [26] . hence invalidating their conclusions. The parity 
violation arises precisely from the presence of the Chern- 
Simons term. 



B. Searching for marginally stable modes 

Having understood the possible source of instabilities 
stemming from the presence of the Chern-Simons term, 
we explain here how we searched for the marginally stable 
modes. In the non-antipodal case the Lagrangian for 
marginally stable modes is still of the form ( 64 ), but now 



the functions a(z) and b(z) are more complicated than 
(65) and depend also on x^z). The Chern-Simons term, 
as it is topological, remains the same. The equations of 
motion following from the Lagrangian ( 64 ) read 



~ h » { z) + ^h'(z) m+ii*kA> (z) ~ = o 

a(z) a{z) 

As explained earlier, the onset of the tachyonic instability 
is usually signaled by a marginally stable mode which is 
normalizable at both asymptotic regions. We search for 
such modes in two different ways. First, we fix the value 



of h(0) by setting it to unity and solve (67) from z = to 
z = oo and from z = to z = - oo for various values of 
h'(0). This approach misses solutions which have h(0) — 
0, so we also need to fix the value of h'(0) (also to unity) 
and solve for the two ranges of z for various values of h(Q) . 
In parity preserving Lagrangians this would correspond 
to looking for even and odd solutions. 



C. Scanning background solutions 

After thorough searches of marginally stable modes for 
various asymptotic separations of D8-brane and various 
baryon charge densities we conclude that in the model 
we study there are no such modes because the electric 
field on the D8-brane is too small. To make the case, 
we plot the values of = h(oo) and = h(— oo) as 
functions of momentum k. A normalizable mode will be 
the one for which and Hr vanish simultaneously. On 
Fig. [6] to [10] we show the results for the antipodal case at 
normalization chosen so that h(0) = 1 for various values 
of the baryon density E and the infrared parameter h'(0). 
Note that scanning for h'(0) > is enough because of the 
symmetry of the equations of motion (flipping the sign of 
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FIG. 6. Asymptotic values of the gauge field perturbation h 
and z — —oo ihh given by the solid black curve) and z = oo 
{Jir given by the dashed gray curve) as functions of the mo- 
mentum of a mode for the baryon density E = 1000, h'(Q) = 
and the antipodal embedding. No marginally stable mode is 
present. The inset plot shows that functions quickly achieve 
large values. Note that the plot is symmetric with respect 
z -H> —5 only because h'(0) = 0. 



h'(0), flipping the sign of k and exchanging z with —z). 
In the inset we show the same plots but for a wider range 
of k. 

We clearly see that the qualitative nature of the plots 
is the same and with the increase of h'(0) just the over- 
all scale of the plots increases. As it is visible on the 
inset plots, hi, and fiR become very large already for in- 
termediate values of momentum k. Thus, if there were 
marginally stable modes they would manifest itself at rel- 
atively small values of k. Although it is possible to make 
the mode normalizable at one of the asymptotic regions 
of the D8-brane, we did not find any example where it is 
normalizable on both ends neither in the antipodal nor 
in the non-antipodal case. 

The absence of this particular marginally stable mode 
turns out to be a subtlety related to the smallness of the 
radial electric field generated by the D4-brane charge dis- 
tributions on the D8-brane. Indeed, if we treat the mean- 
field action in the bottom-up spirit and treat parameters 
appearing there as free, by increasing the baryon mass 
(/?) five-fold, we do see marginally stable modes of the 

As our 



form (61). This is clearly visible on Fig. 
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mean-field description makes various assumptions, it is 
entirely conceivable that microscopic parameters like the 
baryon mass or the assumed DBI form of the kinetic term 
for the gauge field do not enter the mean-field action 
directly but get renormalized in the process of coarse- 
graining. We comment on this possibility in the con- 
clusions. Note also that the absence of this particular 
marginally stable mode is not enough to claim linear sta- 
bility of the system in question. 



FIG. 7. Negative scan for the marginally stable mode at E = 
1000, h'(0) — 1 and the antipodal embedding, notation as on 
Fig. [6] Note that the plot is no longer symmetric with repect 
toz« —z, as h'(0) / 0. 




FIG. 8. Negative scan for the marginally stable mode at E = 
1000, h'(0) — 1.5 and the antipodal embedding, notation as 
on Fig. [6] 




FIG. 9. Negative scan for the marginally stable mode at E = 
1000, h'(G) — 3 and the antipodal embedding, notation as on 
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FIG. 10. Negative scan for the marginally stable mode at 
E = 1000, h'(0) = 50 and the antipodal embedding, notation 
as on Fig. [6] 
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FIG. 11. The appearance of marginally stable mode (h L / R = 
at k ~ 1.2) when baryon's mass is five-fold increased (E = 
100, h'[0] =3 and /3 = 20tt 2 ). 



VI. VECTOR AND AXIAL MESON SPECTRA 
AT LARGE DENSITIES 

A. Mesons at non-zero baryon density 

The key conceptual issue that we are interested in, is 
understanding whether the chiral symmetry can be re- 
stored at low temperatures by finite density effects. In 
the single flavor case that is being considered in this pa- 
per, the analogue of the chiral symmetry is the axial sym- 
metry. One of the ways the chiral symmetry breaking 
would manifest itself in our case is by the parity-odd and 
the parity-even mesons failing to form doublets. Indeed, 
in the spin-one sector that we focus on in the article, the 
axial and vector mesons differ in masses. For the lowest 



two mesons the relative mass difference (29) is as large 
as 40%, which can be seen directly in Fig!~[2 

At non-zero density there are two effects which might 



influence the spectrum of mesons and both require para- 
metrically large densities of order C(A 2 ). For the first 
one, the DBI action is nonlinear and the kinetic term 
for the gauge field perturbations gets modified by the 
presence of the electric field on the D8-brane provided 
the latter is big enough, i.e. 0(A). This exactly corre- 
sponds to baryon densities of order 0(A 2 ). We are not 
aware of any work that addresses this issue for the electric 
field sourced by the four-dimensional holographic baryon 
densities, which is the relevant dual description of order 
0(A 2 ) baryon densities. However, as a similar effect is 
present also in the case of a large enough density of the 
axial charge, it does not seem to be of particular rele- 
vance when looking for possible holographic realizations 
of the Quarkyonic phase. 

The second effect, which to our knowledge has not 
been explored in the literature earlier, are interactions 
of gauge field perturbations corresponding to the axial 
and the vector mesons with the holographic density of 
the baryons. Such a possibility follows from the para- 
metrically large baryon density that we consider. 

In the following, we will analyze the influence of both 
the effects on the masses of the lowest (transverse) vector 
and axial mesons. 



B. The relevance of direct interactions of the 
holographic mesons and baryons 

Let us start with the analysis of how various terms 
in the action scale with N c and A. This will determine 
the parametric regime in the baryon density, in which the 
meson-baryon interactions significantly influence the vec- 
tor and axial meson spectrum. Note that in this subsec- 
tion we will be operating on the original, i.e. not rescaled 
quantities, as we will be discussing the microscopic pic- 
ture. 

The mesons of interest are perturbation of the U(l) 
gauge field on the world volume of the D8-brane. The 
kinetic term for the gauge field, after expanding the DBI 
action to quadratic order, has a factor of 7V C A in fronlp] 
In the vacuum there are no other terms that can com- 
pete with it at quadratic order in the gauge field pertur- 
bations S A and hence the overall scaling is irrelevant for 
the masses of mesons. At non- vanishing baryon density, 
the holographic baryons will couple to the gauge field 
perturbations. Then the factor of (inverse of) iV c A will 
control the strength of electromagnetic interactions with 
single baryons, schematically 



(68) 



9 This can be seen from the equation |[TTJ after Kaluza-Klein re- 
duction on S 4 and expanding it to quadratic order in the field 
strength. 
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where Sj is the current made out of the holographic 
baryons. This equation is nothing but an approximate 
and simplified version of (51). As the baryon charge 



scales as -/V c , the current Sj schematically reads 
Sj = N c nSu. 



(69) 



no restoring force generated from interactions with other 
charges from the lattice. We can estimate the magnitude 
of the restoring force by expanding the interaction energy 
of a single baryon with others (notation the same as in 
the equation (43 1) up to quadratic order in its displace- 
ment Sx 



We assume that perturbations of interest do not affect 
the number density n, which will indeed turn out to be 
the case. 

As meson perturbations are, in the single flavor case, 
the electromagnetic waves on the D8-brane, there will be 
a Lorentz force acting on the holographic baryons. We 
can schematically write Newton's second law ma = qE 
as 



N r XdSu = NrdSA, 



(70) 



where iV c A on the LHS comes from the baryon mass, 
whereas N c and dSA come respectively from the baryon 
charge and the electric component of the electromagnetic 
wave on the D8-brane. Although, as follows from ( [70| , 
baryon acceleration^ and hence baryon velocity!""! will 
be 1 / A-supressed 



<y« = o(i/A), 



(71) 



for C(A 2 ) number densities n the current Sj will con- 
tribute to the equation for meson perturbation ( 68 ) at 



the same order of large iV c and large A scaling as the 
kinetic term. This implies that the electromagnetic re- 
sponse of baryons may significantly affect the propaga- 
tion of electromagnetic waves on D8-brane and hence the 
meson spectrurrp] 

One can clearly see that in the case of the marginally 
stable modes studied in previous section, there will be no 
Lorentz force acting on baryons. This is because baryons 
are at rest in equilibrium and the marginally stable mode 
(61 1 does not generate an electric field. This justifies ne- 



glecting the dynamics of charges in the analysis of the 
previous section. When it comes to the meson perturba- 
tions, these necessarily introduce an electric field in one 
of the spacelike directions. Such an electric field will ac- 
celerate charges according to (70). However, once one of 



(restoring force) • Sx 



E 
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A d? 
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(72) 



leading to 



restoring force - 



Nc_ 
d 4 \ 



Sx fa —N r X Sx. 



(73) 



This implies that the equation ( 70 ) we started with needs 
to be supplemented with the additional term (73). Tak- 
ing into account that 



Su = Sx, we end up with the 
schematic expression resembling a forced harmonic os- 
cillator 



N c X(Sx + luqSx) = N C SA. 



(74) 



Note that for velocities 0(1/A) (and hence for displace- 
ments of the same order), both sides of (74) are of the 



same order in A, which means that the lattice structure 
does not prevent the holographic baryons from being in- 
fluenced by electromagnetic waves representing the holo- 
graphic mesons. As all factors of A cancel both in (74) 



and (68), we conclude that interactions of the baryonic 



medium with the mesons will influence the properties of 
the latter. 

At this moment we do not have a precise form of the 
equation ( 74 ) , not to mention the precise value of the 
coefficient loq. In particular, note that if Sx denoting 
displacement in one of the spacelike field theory direc- 
tions depends only on time, then there is no restoring 
force as the lattice as a whole moves uniformly. This 



implies that Sx must enter restoring force formula ( 73 ) 



the charges is displaced, one may wonder whether there is with the derivative in the radial direction, which is the 



reason why ( 74 ) only resembles a forced harmonic oscil 
lator. Note also that such a term might appear in the 
effective description (48) as a derivative correction. 



10 If there is an electric field associated with perturbation. 

11 Assuming oscillations with O(l) frequency, as is the case with 
meson masses in vacuum. 

12 One can easily see that the form of the equations ( |68[ t — ( f^0| ) is 
not affected by simultaneous rescalings of 5 A and Su by a factor 
of A. This is the reason why meson-baryon interactions of the 
form j68t-j70| are captured in a precise fashion by the rescaled 
action~ |l8[ l + \48\ . Although there we were arguing that we 
shall never be interested in studying macroscopic velocities in 
spacelikc directions, because of the scaling argument elucidated 
above, the action | |18| l + ( |48| can be nevertheless used to (partly) 
calculate the equations of motion for meson perturbations in the 
presence of interactions with baryons. 



In the following, we will explore the implications of 
interactions of the mesons with the baryons completely 
neglecting the restoring force. This could be justified if it 
turned out that ujq in the cases of interest was much lower 
than the mass of lowest vector or axial meson. We do not 
know this at the moment. Our motivation is rather to 
show that density-enhanced holographic meson-baryon 
interactions may lead to very interesting physical effects, 
which provides a very strong motivation for more detailed 
future explorations of the phase diagram of the Sakai- 
Sugimoto model at 0(A 2 ) baryon density. 
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C. Gauge field perturbations describing vector and 
axial mesons at 0(X 2 ) baryon density 

In the following, we will search for the normalizable 
gauge field perturbations on the D8-brane of the form 
(26). As explained in subsection II A the ansatz (26) 



corresponds to zero momentum axial and vector meson 
currents in the x 3 direction. Their frequencies have direct 
interpretation in terms of meson masses. In the vacuum, 
the differences in the masses of the corresponding vector 
and axial mesons come from the broken axial symme- 
try 



Here we will reexamine perturbations ( 26 1 , but now 



in the presence of a parametrically large baryon density, 
as considered in section |II A| The goal is to obtain the 
frequencies of the normalizable modes, which we under- 
stand as masses of mesons in the dense baryonic medium, 
and compare them with the masses of the corresponding 
mesons in the vacuum. If the chiral (axial) symmetry is 
(approximately) restored, then the masses of the corre- 
sponding axial and vector mesons will be (approximately) 
equal. 

The notion of chiral symmetry restoration that our de- 
scription possibly allows for is the one in which pertur- 
bations on the _D8-brane and D8-brane (approximately) 
decouple because of the interactions with the holographic 
baryons. 



Indeed, for perturbations of the form (26 1 this might be 



in principle possible, as the ansatz ( 26 ) leads to a time- 
dependent electric field in i 3 direction. The electric field 
accelerates baryons leading to a time-dependent, but uni- 
form, current in the i 3 direction. This current backreacts 
on the electric field through the standard electromagnetic 
coupling present in the action (481, which (intuitively) 



lowers the original electric field in complete analogy with 
the London equation in case of an electromagnetic field 



13 A very interesting question that we are not going to analyze in 
this article is the fate of the holographic mode corresponding 
to the pion (rf meson) at parametrically large baryon densities. 
If the pion remains massless, at zero momentum its frequency 
should vanish. The pion corresponds to a perturbation of the 
form A M = i/>(z)9 ) j(/)(x), where ip(z) is antisymmetric and ap- 
proaches non-zero and opposite constant values on the bound- 



aries [15] . At zero momentum, i.e. (p(x) 



such a mode 



generate an oscillating radial electric field and, as our system has 
broken translational invariance in the radial direction, studies of 
this mode in the presence of the dense baryonic medium seem to 
be more difficult than that of the modes corresponding to vector 
and axial mesons at zero momentum. If there is a non-zero fre- 
quency for which this mode exists, than pion effectively becomes 
massive, which would suggest that the chiral (axial) symmetry 
might be restored. It would be very interesting to analyze this is- 
sue in more detail, but the following argument suggests the pion 
will remain massless: if the spectrum does not jump discontinu- 
ously as we increase the density, the pion will still correspond to 
a perturbation of the form = i/)(5)S M 0(x). This perturbation 
does not change under <j>(x) — > <j)(x) + const, and therefore this 
must be a symmetry of the system, which in turn can only be 
the case if the pion is massless. 



trying to penetrate the interior of the superconductor. 
The full picture will turn out to be more intricate be- 
cause of the z-dependent background geometry and the 
charge density. 

Of course, in contrast to the black hole embedding 
when the chiral symmetry is exactly restored, for mesons 
of sufficiently large mass the holographic baryon density 
will not be very relevant. We will revisit this issue in 
the following subsection by doing a Schrodinger equation 
analysis of meson perturbations searching for the poten- 
tial barrier emerging from bulk meson-baryon interac- 
tions and in this subsection we will present numerical 
results of the lowest part of the vector and axial meson 
spectrum as a function of the asymptotic separation of 
D8- and D8-branes and the baryon density B. 



The equations of motion for the perturbation ( 26 ) will 



follow from taking the functional derivative of the action 
(18) + (48) with respect to 60,3 (t,2) and 8v,3(t,z). After 



going to frequency space, the equation for 80,3(2) takes 
the schematic form 



5~4(z) 



Ci (z)6a' 3 (2) + Q 2 C^ (2)8a s (2) 



w 2 ' 



C^ w) (z)Su 3 (z)=0, 



(75) 

where the exact form of Cq(2), Cq(2) and C\(2) is not 
very illuminating and can be found in the appendijj^] 
Note also the presence of the term w 2 Cq(2) which was 
absent in the vacuum case and here encodes the interac- 
tions of mesons with baryons. This term will lead to new 
physical effects. 



The equation for 6113(5) takes the form of (50dl 



r yw(2) 2 60.3(2;) = — 2X(z) Sv,3(z) 



with A given by 



\ = -p~e , \i + tf/ i2 w 2 . 

2 V 



(76) 



(77) 



The formula above is the curved space analogue of the 
London equation and leads to a decoupled equation for 
the meson perturbation 

65,1(2) + C 1 (2)5a' 3 (2) + cj 2 C*^ ) (2)8a 3 (z) + 

C (2)8aa(2) = 0. (78) 

The form of the coefficient Cq (2) can be found in the 

appendix and together with C\(2) and Cq W \2) it encodes 
the form of the embedding, the background electric field 
and the interactions of mesons with baryons. It can be 
easily verified that the equation for meson perturbations 
is symmetric with respect tozf> — 2 and hence it is suffi- 
cient to solve it for 2 between and 00. Imposing Dirich- 
let boundary conditions at infinity leads to the spectrum 
of mesons at 0(X 2 ) densities. 



14 Actually, we skip the form of Cq 1 "'^), as in the master equation 
( |78| it appears in a combination with other term and in the 
appendix we only provide the form of the resulting expression 
denoted by Cq(z)- 



20 



D. Meson spectra and the effective potential for 
meson perturbations 

The presence of density-enhanced direct holographic 
meson-baryon interactions changes the form of the equa- 
tion for meson perturbations, hence influencing the spec- 
trum of the axial and vector mesons. In order to under- 
stand better the features introduced by the interactions 
in question, it will be extremely useful to rewrite the 



equation for perturbations ( 78 1 in the Schrodinger form, 
i.e. 



2(fcr2 



(79) 



where the energy E of the wave function is the square 
of the rescaled frequency U) and hence is proportional to 
the square of meson mass m (27 1 



E = = ~™ 2 /M 2 KK , 



(80) 



vE^cr) is related to 03(f) and a is a new radial coordinate 
that leads to the canonical kinetic term in ( 79 ) . Simple 



manipulations on ( 78 1 allow us to obtain the potential 
term V. For simplicity V can be expressed as a function 
of z instead of a 



V S chr{z) 



Co(i) Ci(S) 2 5C^'(zf 



C^\z) 4Cr(S) 



C[(z) 



4C^ ] (z) 2 '' 



(81) 



where Co (5), C (z) and C\(z) are as (78). 

Because of the way the background spacetime is 
warped, V needs to have the form of an infinite potential 
well. This feature gives rise to an infinite tower of bound 
states, which are simply the vector and axial mesons. 
However, the presence of a non-trivial embedding, non- 
zero background electric field and most importantly holo- 
graphic interactions of the baryons with the mesons may 
lead to new features of the potential having a direct im- 
pact on the spectrum of axial and vector mesons. 

In the following we will analyze in detail two cases, 
which based on more exhaustive studies, we think are 
representative ones. 

Fig. [12] shows the Schrodinger potential in vacuum 
(upper plot) and at very large density (E = 1000, lower 
plot) for the antipodal embedding. The green dashed 
and red dotted lines denote the squares of vector and 
axial meson masses for the three lowest mesons. As an- 
ticipated earlier, the potential has the structure of a well 
supporting infinitely many bound states. In the vacuum 
case, the potential is simple and does not have any finer 
structure than the master infinite potential well coming 
from the warping of the target spacetime. 

If one looks at the situation at non-zero density, one 
sees that the potential develops substructure in the form 
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FIG. 12. The Schrodinger potential and the three lowest ax- 
ial and vector mesons in the antipodal case (L/St — 0.5) 
at vanishing (top plot) and large (bottom plot, E = 1000) 
density. Green dashed and red dotted lines correspond to 
the squares of masses of the lowest three vector and axial 
mesons respectively. The most interesting finding is the ap- 
pearance of an emergent potential well due to the holographic 
baryons on the D8-brane. It can be clearly seen, that in the 
presence of a dense baryonic medium, the lowest axial and 
vector mesons are bound states not of the master potential 
well, which is present also in the vacuum, but rather of the 
emergent one. This leads to a tempting connection with a 
quasiparticle picture in which some of the lowest in-medium 
mesons are collective excitations of underlying Fermi surface. 



of another, finite potential well, which is entirely due to 
the presence of charges on the flavor brane. The immedi- 
ate observation is that the lowest vector and axial vector 
mesons arise then as bound states of the emergent po- 
tential well, rather than the original one. This strongly 
suggests the analogy with the quasiparticle picture in 
which low energy long-lived excitations, when expressed 
in terms of the underlying microscopic description, turn 
out to be collective modes of a system. Note also that 
the excitations are automatically color singlets, as is the 
case in the Quarkyonic phase. 

The same feature, that the lowest vector and axial 
mesons arising as bound states of the emergent poten- 
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FIG. 13. The Schrodinger potential and the three lowest ax- 
ial and vector mesons for a very non-antipodal embedding 
(L/St — 0.05) zero (top plot), medium (middle plot, E = 100) 
and high (bottom plot, E = 1000) density. Because of the 
small potential barrier residing in the vicinity of z = 0, the 
emergent potential well does not give rise to bound states 
and effectively acts as a potential barrier. On the other hand, 
the emergent barrier leads to an approximate chiral symmetry 
restoration for the two lowest states, as their relative mass dif- 
ference given by ( |29[ ) reduces from about 40% in the vacuum 
reduces to 10% in this case. 




FIG. 14. The Schrodinger potential in the antipodal case in 
the vacuum (solid curve) and at densities E — 100 (dashed 
curve), E = 500 (dash-dotted curve) and E = 1000 (dotted 
curve). One clearly sees that the emergent potential barrier 
arising because of the charges on the D8-brane will affect more 
and more of the lowest lying axial and vector mesons as the 
density increases. 



tial associated with the presence of the four-dimensional 
density of holographic charge can be seen already at low 
rescaled densities E. This becomes apparent from Fig. 
|14| showing the Schrodinger potential in the antipodal 
case for the vacuum and three representative charge den- 
sities and Fig. [15] showing the masses of the the three 
lowest axial and vector mesons as functions of density, 
also in the antipodal embedding. 

We are not aware of any similar results in the gauge- 
gravity duality, but we also want to stress at this point 
that our studies are preliminary and clearly, as explained 
in the beginning of this section, do not take into account 
all in-medium phenomena. 

Setting this aside, we also want to point out another 
possibility that the density-enhanced meson-baryon in- 
teractions offer, that the emergent potential well from 
the point of view of the master potential well might also 
act as a potential barrier. This clearly does not happen 
in the antipodal case, where the emergent potential well 
supports a few bound states, but the asymptotic separa- 
tion of the D8- and D8-brane L is a free parameter and 
we can investigate whether the situation changes when 
we set it to some different value. 

It turns out, that in the non-antipodal case in the vac- 
uum there is already a small potential barrier in the mid- 
dle of the master potential well, as seen on Fig. 13 (top). 
Again, the charge on the D8 will result in the creation of 



an emergent potential well, but now this small potential 
barrier will not allow for any bound states inside it. It 
is thus a genuine potential barrier which offers the pos- 
sibility of forming chiral doublets - (almost) degenerate 
states residing on the two sides of the barrier. Indeed, 
looking at the two lowest states (higher excitations have 
energies larger than the height of the barrier) we see that 
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FIG. 15. Masses of the three lowest vector (top) and axial 
mesons (bottom) as a function of rescaled baryon density E 
for the antipodal embedding of D8-brane. One can see that 
within the adopted approximations the masses obtain signifi- 
cant contributions from interactions with the dense baryonic 
medium. Saturation of masses with density follows from the 
fact that the charge distribution is not altered by adding new 
layers, which leads to the Schrodinger potential as on Fig. |14| 



the relative ratio of masses of axial and vector mesons 
(29), which, as shown on Fig. [2j in the vacuum is of 
order of 40%, in this case can become as small as 10%. 



We interpret this result as approximate chiral symmetry 
restoration in the lowest part of the spectrum. Clearly, 
for higher excitations the potential barrier plays a less 
important role. 

The last point of our analysis is to check what is the 
precise role of the density-enhanced direct holographic 
meson-baryon interactions on the form of the potential. 
Fig. [16] shows the Schrodinger potential with this ef- 
fect turned off by hand and one can clearly see that the 
most interesting feature present before, i.e. an emergent 
relatively deep potential well, is absent now. Although 
our analysis is far from complete, based on our findings 
we think that capturing, the so far neglected, density- 
enhanced direct holographic meson-baryon interactions is 
crucial for getting the physics of dense holographic QCD 
models right. 

To conclude this section, we would also like to stress 
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FIG. 16. Schrodinger potential in the antipodal case for den- 
sity E = 1000 in the absence of direct holographic meson- 
baryon interactions, i.e. Co (5) from (78 1 encoding the holo- 



graphic baryons' response to the electric field perturbations 
is set, by hand, to zero. Comparing with the second plot 
within Fig. [12] we see that the emergent potential well gets 
much more shallow and lacks the barriers present when di- 
rect holographic baryon-meson interactions are taken into ac- 
count. The discontinuity of the potential at the boundary of 
charge distribution follows from the discontinuous first deriva- 
tive of background radial electric field. 



that in all examples we analyzed we never found states 
that we could not connect continuously to the ones en- 
countered in the vacuum, as is neatly depicted in Fig. 
(15) for the three lowest vector and axial mesons in the 
This means that the bound states we 



antipodal case, 
found were never genuinely new mesons, but rather the 
ones already existing in the vacuum with properties sig- 
nificantly altered by the presence of the dense baryonic 
medium. If the effect we describe persists in more de- 
tailed microscopic studies, than it might have interesting 
consequences that we discuss in the last section. 



VII. SUMMARY AND CONCLUSIONS 

A. General motivation 

In this paper we have considered the single flavor Sakai- 
Sugimoto model at zero temperature and finite baryon 
density. Our main motivation was to look for evidence 
of the existence of the so-called Quarkyonic phase [5] . In 
the Quarkyonic phase, baryons are so tightly packed that 
quarks no longer belong to a particular baryon, rather, 
one effectively has a quark Fermi surface. Despite that, 
the excitations of the systems are, as in the vacuum, color 
singlets, which motivates the name: Quark (from having 
the quark Fermi surface) - yonic (from excitations be- 
ing colorless particles including baryons). The presence 
of the quark Fermi surface has profound consequences 
for the properties of the Quarkyonic phase: chiral sym- 
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metry is expected to be restored [5] and the system is 
unstable against the formation of Quarkyonic chiral spi- 
rals [§] [TU] . The latter break the chiral symmetry again, 
together with the translational symmetry, but the chi- 
ral condensate oscillates in space around zero with the 
wavelength related to the chemical potential. 

Since discussions of the Quarkyonic phase in the QCD 
literature are phenomenological in nature, it is interest- 
ing to see to what extent the gauge-gravity duality can 
shed light on the existence and properties of this phase. 
This was the chief motivation for our studies and be- 
low we summarize our results detailing once more the 
assumptions made and listing the most promising open 
directions. 

We began our analysis by motivating in section [TTT| why 
low temperatures and baryon densities of order 0(X 2 ) 
- represented holographically by four-dimensional lat- 
tices of the holographic baryons having extent both in 
the field theory directions and the holographic direction 
[3"T] - might correspond to the Quarkyonic phase in the 
holographic model of interest. We listed three reasons 
that are characteristic of this regime of parameters: the 
presence of strong electromagnetic interactions between 
holographically represented baryons, which might be in- 
terpreted as a dual manifestation of the formation of a 
quark Fermi surface [31]; the possibility of spontaneous 
breaking of translational invariance because of instabili- 
ties triggered by the Chern-Simons term [5S1IM]; density- 
enhanced interactions of the holographic mesons with 
baryons, which significantly modify the meson spectrum. 
Although the first two effects were already anticipated 
in the holographic literature, albeit never directly in the 
context of the Quarkyonic phase, the third one, to the 
best of our knowledge, is a completely new phenomenon. 



B. Coarse-grained description of holographic 
baryonic medium 

The chief difficulty that studies of very dense holo- 
graphic QCD are faced with lies in constructing holo- 
graphic arrays of baryons corresponding to an C(A 2 ) 
baryon density in the dual field theory. A fully micro- 
scopic solution of the problem would amount to mini- 
mizing the energy of a large number of point-like charged 
particles interacting through DB1 electromagnetism, put 
in a finite volume in the field theory directions and in 
an external gravitational field associated with the holo- 
graphic direction. Although the first attempts to con- 
struct such states were already taken in (35], the solu- 
tions obtained there consisted of only a few layers of the 
holographic baryons in the holographic direction, which 
is superficially too small a number to observe any of the 
effects that we mentioned. 

Because of these difficulties, following [31], we adopted 
a simple coarse-grained description in which we minimize 
the energy generated by the continuous distribution of 
the holographic baryons. This way of describing the sys- 



tem turned out to be very similar to the electron star 
construction of [36] and we drew this analogy in our con- 
siderations, making contact between the two approaches 
in section |IV[ The effective Lagrangian originally pro- 
posed in [31 1 and used by us is very simple and postu- 
lates a DBI form of the kinetic term for the mean-field 
U(l) gauge field and the D8-brane mean-field embedding 
function and a minimal coupling of the density of the 
holographic baryons to the mean-field U(l) gauge field. 
In the model we also assume that the energy density of 
the holographic medium is that of a free system. 

All these assumptions might in principle turn out to be 
violated in a true top-down model of the very dense holo- 
graphic baryonic medium. In particular, the approach we 
followed does not take into account oscillations of the mi- 
croscopic gauge field and the embedding function over the 
lattice spacing scale. It is entirely conceivable that such 
effects will modify the form of the effective Lagrangian, 
e.g. modify the equation of state or the DBI form of 
the kinetic term for the mean-field U(l) gauge field. Al- 
though we expect the qualitative, rough features of our 
findings to hold also for truly microscopic charge distri- 
butions, quantitative answers might be significantly dis- 
torted by the simplifications made. This might partly 
justify promoting the microscopic quantities appearing 
in the coarse-grained model, like the mass of the holo- 
graphic baryon and the parameters ft and 7 introduced in 
section [IV B[ to free parameters in the bottom-up spirit. 

C. Antipodal and non-antipodal embeddings 

The original Sakai-Sugimoto model with two D8- 
branes located antipodally in the circular direction 
has a generalization with the asymptotic separation of 
D8-branes being an extra parameter having no direct 
analogue in QCD. The phenomenology of the Sakai- 
Sugimoto model does depend on this parameter, denoted 
here by L. For example, the phase diagram is affected 
by L |50j , and Yukawa interactions between baryons and 
scalar mesons [51] depend on L and will vanish in the 
antipodal case due to the extra Zi symmetry. 

This is the reason why the first new result delivered in 
our studies was generalizing the coarse-grained descrip- 
tion of [3T] to non-antipodal embeddings. Interestingly, it 
turned out that within the approximations made by us, 
the properties of vector and axial mesons in the dense 
baryonic medium depend rather strongly on L, which a 
posteriori is another motivation for this modification of 
the original Sakai-Sugimoto model. 



D. Holographic charge distribution and the quark 
Fermi surface 

Our coarse-grained description naturally leads to four- 
dimensional densities of the holographic baryons , as seen 
also in the more microscopic calculations of [32]. This is 
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the most distinctive feature of the considered model at 
C(A 2 ) baryon densities, as opposed to lower densities, 
and, together with the common lore of identifying the 
holographic direction with the energy scale in the dual 
field theory, one sees that the holographic charge distri- 
bution is somewhat reminiscent of a Fermi sea. In fact, it 
is tempting to identify the edge of the holographic baryon 
distribution directly with a Fermi surface, as was done 
by the authors of [31] . Some words in caution are in or- 
der here. In the first place, although it is true that in 
the ground state the lowest lying states that are avail- 
able to the baryons are being occupied, it is important 
to remember that it is not the Pauli exclusion principle 
but the Coulomb repulsion which prevents baryons from 
occupying the same state in the bulk. Moreover, the 
states are being occupied in the coordinate (gravity dual) 
space rather than momentum space, and the baryons are 
strongly interacting so a naive weakly coupled second 
quantized point of view is inappropriate. Having said all 
that, we think that the best indicator of the presence of 
a quark Fermi surface is the very strong overlap between 
the baryons in the dual field theory and at the same time 
the strong interactions between their holographic coun- 
terparts. Although the first feature is present at much 
lower, even 0(1) baryon densities, the latter arises only 
for a parametrically dense medium. 



E. Chern-Simons term-induced instabilities 

Apart from possessing a quark Fermi surface, one of 
the characteristic features of the Quarkyonic phase is the 
appearance of unstable modes breaking translational in- 
variance. In the holographic literature it is known that 
the presence of the Chern-Simons coupling for electro- 
magnetic/DBI interactions together with the radial elec- 
tric field may lead to the appearance of unstable modes 
breaking translational invariance [33-35] . So far this 
mechanism, which appears in a couple of variants [541 142 j . 
seems to be the only one known within the gauge-gravity 
duality that might lead to time-independent modulations 
in holographic media. One of the indications of the 
instability is the presence of marginally stable modes, 
i.e., time-independent normalizable perturbations carry- 
ing non-zero momentum. A mode of that type, involving 
a helical modulation in baryon or axial currents, is known 
to appear in the Sakai-Sugimoto model in the deconfined 
phase [29] , as well as at zero temperature but large axial 
chemical potential [30] . 

As the appearance of this particular mode in the Sakai- 
Sugimoto model superficially seemed to be rather robust, 
we performed a thorough search for it for various asymp- 
totic separations of the flavor branes, as well as various 
charge densities. Details, including interesting technical- 
ities, can be found in section [V[ It turned out that within 
the adopted coarse-grained description we could not find 
such a mode, which we associated with the smallness 
of the radial electric field generated by our distributions 



of the holographic baryons. However, upon promoting 
the microscopic baryon mass appearing in our coarse- 
grained lagrangian to a free parameter, we started seeing 
marginally stable modes for masses about five-times big- 
ger than the microscopic ones, as can be seen in Fig. 11 
The chief difference between the two models is that the 
radial electric field generated by the charges on the fla- 
vor brane gets a couple of times bigger in the latter case. 
This might be regarded as an indication that our coarse- 
grained model needs to be modified in order to capture 
the physics of interest. It is also important to stress that 
while searching for instabilities we did not consider the 
most general ansatz and, hence, there might be unstable 
modes that evaded our analysis. We leave this issue as 
an open problem. 

Our instability analysis in section |V| referred to insta- 
bilities in the gauge field, whose time-like component is 
dual to the left and right quark number operators ip^L 
and ipt^R. Potential instabilities of this type would 
spontaneously break translational invariance but not af- 
fect chiral symmetry. The Quarkyonic chiral spiral in- 
stability on the other hand, refers to spatially dependent 
expectation values for both the chiral condensate as well 
as operators of the form tploli"^ which also mix left and 
right quarks. It is possible that there exists a link be- 
tween the two instabilities, and it would be interesting 
to explore this further. It might well be that evaluating 
the chiral condensate in the presence of the condensed 
Cherm-Simons term-induced instability gives an expres- 
sion similar to the one appearing in the Quarkyonic chiral 
spiral. 

It is also worth pointing out here that if there is chi- 
ral symmetry restoration in the Quarkyonic phase, the 
presence of the Quarkyonic chiral spiral would break the 
chiral symmetry again but with the condensate averag- 
ing to zero over sufficiently large scale. To reproduce 
the phenomenology of the Quarkyonic phase, we should 
therefore either find an unstable phase with chiral sym- 
metry restoration, or a stable phase with broken trans- 
lational invariance and broken chiral symmetry with the 
chiral condensate oscillating in space around zero. 



F. The spectrum of mesons at parametrically large 
baryon densities 

As we have emphasized several times, it is only at 
baryon densities of order C(A 2 ) that interactions of the 
flavor brane gauge field perturbations with the holo- 
graphically represented baryons can significantly influ- 
ence the spectrum of the normalizable modes of the gauge 
field and hence the masses of the corresponding vec- 
tor and axial mesons. This happens because the holo- 
graphic baryons are charged particles which accelerate 
in the presence of the electromagnetic field perturbation 
and in this way lead to a non-trivial dynamical electric 
current coupled to DBI electromagnetism on the flavor 
brane world volume. 
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Naively such an effect should not be visible, as baryons 
are heavier than their charge by a factor of A and their 
acceleration due to the electric field associated with 
the gauge field perturbation and hence their velocities 
and displacements from equilibrium positions are mass- 
suppressed. However, the large density of the holographic 
baryons compensates the smallness of their velocities and 
- this is the crucial observation here - the current pro- 
duced by the holographic baryons yields a contribution 
to the equation of motion for the gauge field perturbation 
which is as important for the dynamics as the gauge field 
kinetic term. This effect is analogous to the one described 
by the London equation for electromagnetic waves try- 
ing to penetrate a superconductor and can be under- 
stood by the effective mass the gauge field acquires inside 
the medium. This analogy superficially offers a natural 
mechanism for the chiral (axial) symmetry restoration - 
a massive gauge field might not be able to penetrate the 
charge distribution leading to an effective decoupling of 
modes residing on both sides of it. 

Although this is what we originally expected to find, 
the full story turned out to be more complicated. The 
modifications arise in the first place because the gauge 
field perturbations propagate in a curved space-time and 
try to penetrate an inhomogeneous charge distribution 
and, secondly, the holographic baryons form a lattice and 
hence one might expect the presence of a restoring force 
acting on a baryon outside its equilibrium position. The 
simple analysis in section (VI) gave an estimate of the 



magnitude of the latter effect, and it turned out to con- 
tribute at leading order, at the same order as the kinetic 
term for the gauge field and as the contribution from 
the current of the holographic baryons . Understanding 
the consequences of this potentially very important effect 
requires the knowledge of the details of the holographic 
baryon distribution that we clearly lack and in our anal- 
ysis we completely neglected it. Our motivation was to 
show the existence of new phenomena at C(A 2 ) baryon 
densities and to motivate a further study, rather than to 
incorporate all the intricate details right from the outset. 

A neat way of understanding the eigenvalue prob- 
lems of the type encountered here is to write the master 
equation for perturbations in the form of a Schrodinger 
equation. In this picture the potential appearing in the 
Schrodinger equation allows one to understand the qual- 
itative features of the spectrum obtained through direct 
calculations. 

In the vacuum for the antipodal embedding the poten- 
tial has the form of a single well that is infinite and sup- 
ports infinitely many bound states (top plot on Fig. 12 ). 



In the non-antipodal case a small substructure emerges, 
i.e., a small potential barrier whose height however does 
not affect significantly the properties of the mesons (top 
plot on Fig. 13 ) . In particular, for a wide range of asymp- 



of chiral symmetry breaking in large- N c gauge theories. 
Put another way, if the chiral symmetry were restored, 
then the vector and axial mesons would appear in dou- 
blets, and in particular would have had equal masses. 



The spectrum of vector and axial mesons turn out to 
change significantly if we go to 0(X 2 ) densities and in- 
clude the effects coming from the curved space analogue 
of the London equation. The most striking feature that 
arises is the appearance of an emergent potential well 
inside the master potential well. The lowest vector and 
axial mesons appear then as the bound states of the emer- 
gent potential rather than the one existing already in the 
vacuum. This can be very clearly seen in Fig. |12| A 
phenomenon of that type, to the best of our knowledge, 
has not been yet encountered in the studies of any of the 
existing holographic models. 

The Schrodinger potential at 0(A 2 ) gets even more in- 
teresting when we consider non-antipodal embeddings. 
We already discussed that the emergent potential barrier 
can support bound states by itself and that in the non- 
antipodal embeddings the geometry leads to the appear- 
ance of a small potential barrier in the Schrodinger po- 
tential in the vacuum. At large enough density these two 
features combine leading to an effective potential barrier 
inside the master potential well. This potential barrier 
screens the wave functions on both sides and approxi- 
mately restores the chiral (axial) symmetry between the 
lowest vector and axial mesons as can be seen in Fig. [13] 
However, higher states have energies well above the emer- 
gent potential barrier and the chiral symmetry cannot be 
restored for them by means of this mechanism. 



It is important to stress that we found no genuinely 
new states as a result of our analysis of dense holographic 
QCD - the axial and vector mesons at C(A 2 ) baryon den- 
sities could be continuously connected to the ones ap- 
pearing already in the vacuum as visible in Fig. 
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totic separations the relative ratio of the masses of the 
lowest vector and axial vector mesons turns out to be 
about 40% (Fig. [2]). This can be attributed to the axial 
symmetry breaking, which is the single flavor analogue 



Concluding this point, we found evidence for signifi- 
cant modifications of vector and axial meson spectra at 
baryon densities of order 0(\ 2 ). Our analysis was clearly 
incomplete as we neglected the restoring force acting on 
baryons within the holographic lattice. This might be 
justified if it turned out that the numerical coefficient 
associated with this effect was an order of magnitude 
smaller than the one driving the coupling of the baryon 
current to the gauge field, but we have no evidence for 
this. For the moment we do not have any further detailed 
insight to contribute to this issue. However, the message 
we want to convey is that there is a set of new phenomena 
that affect the spectrum of mesons at baryon densities of 
order C(A 2 ) that seem to have been overlooked in the lit- 
erature on this subject, but their precise effect still needs 
to be worked out in detail. 
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G. Chiral condensate and chiral symmetry 
restoration 

A standard way to examine whether the chiral sym- 
metry is broken or not is examine whether the chiral 
condensate ("ipip) — (tp^tpE + tPriPl) vanishes or not. In 
the Sakai-Sugimoto model, the left and right chirality 
quarks sit at different points on the boundary circle, and 
the gauge-invariant operator which generalizes the chi- 
ral condensate is {ip\P exp(J L A t (1t)^r + L o R). To 
compute its expectation value, one needs to compute the 
regularized area of a world-sheet disc which ends on the 
D8-D8 branes and on the boundary [55] . In the presence 
of a finite density of baryons, and in the mean- field de- 
scription that we have been employing, the computation 
needs to be somewhat modified. The non-trivial gauge 
field on the D8-D8 branes couples to the world-sheet and 
will give rise to an extra contribution 

chiral condensate cx Pexp (^J dzA^j x area(disc). 

(82) 

It might appear that the extra contribution vanishes since 
we chose the gauge Az = 0. However, we should evalu- 
ate the extra Wilson line in the z-direction in the gauge 
in which all gauge fields vanish at infinity. One can then 
see that the extra Wilson line is related to the pure-gauge 
behavior of A^ at infinity in the Ag = gauge. As briefly 
mentioned in section |II A[ and discussed in more detail 
in [T3] , the pure gauge behavior of A^ encodes the expec- 
tation value of the pion field, and we can more precisely 
write 

chiral condensate cx exp(2II(a;' i )// w ) x area(disc), (83) 

where II (x M ) is the pion field and f v is the pion decay 
constant. We therefore see that there is a direct rela- 
tion between the chiral condensate and the expectation 
value of the pion field in our model. In the mean field 
description the pion field vanishes and therefore the chi- 
ral condensate does not vanish. In fact, in the antipodal 
case the geometry of the embedding does not change as 
we increase the baryon density, and therefore the chiral 
condensate does not change either. In the non-antipodal 
case the embedding does change, and so does the chiral 
condensate, but not in a drastic way and so chiral sym- 
metry appears to remain broken. We will return to this 
point below. 

We saw that chiral symmetry became less severely bro- 
ken in the presence of finite density but did not quite 
get restored. It is tempting to speculate about possible 
qualitative mechanisms that could cause effective chiral 
symmetry restoration, as this was one of the main mo- 
tivations for this work. So far, the only mechanism for 
chiral symmetry restoration in these classes of models al- 
ways involves finite temperature physics and a bulk black 
hole. If the preferred phase involves branes extending all 
the way to the black hole horizon, then chiral symmetry 



is restored, otherwise it appears to be broken. So what 
possible mechanisms could one imagine? 

One possibility is that in the presence of the finite 
density medium, gauge fields become effectively mas- 
sive. This effect in our case led to a very intricate form 
of the Schrodinger potential, including in some cases 
the appearance of a potential barrier associated with 
the presence of the holographic baryons, but it never 
caused the gauge field to decay rapidly inside the holo- 
graphic medium. If this would have been the case, then 
the "Fermi surface" would effectively impose Dirichlet 
boundary conditions on the gauge field. Consequently 
the gauge fields on the two sides of the medium would 
become decoupled from each other and chiral symmetry 
would be restored. 

A second possibility that exists is that the backreaction 
becomes important. If the finite density medium would 
cause a large backreaction on the bulk geometry causing 
the geometry to extend further in the IR, then the D8-D8 
system could also extend further into the bulk, causing 
approximate but not exact decoupling between the two 
ends. However in our setup backreaction does not seem 
to be important since it is suppressed by Nf/N c and will 
become important only at 0(N C ) baryon densities. 

Another possibility is that chiral symmetry is not ex- 
actly restored but only approximately on long distances. 
Such a scenario has been proposed in the context of the 
Skyrme model, where one models the dense nuclear mat- 
ter using a Skyrme crystal (see e.g. [53J [M])- In such 
systems, the Skyrme field exp(2n(x A ')// 7 r) which also ap- 
peared in [23 fluctuates rapidly and averages to zero over 
longer distances signaling chiral symmetry restoration. 
To see similar behavior in our case, we would really need 
to take the crystal nature of the ground state into ac- 
count. The closest analogue of the Skyrme lattice in our 
model is a three-dimensional lattice of baryons all located 
at the tip of the D8-brane. As soon as the lattice becomes 
four-dimensional, the connection to the Skyrme lattice is 
less clear. Nevertheless, it would be interesting to investi- 
gate whether in our model a more detailed description of 
the crystal does give rise to a rapidly fluctuating Skyrme 
field exp(2n(x M )// 7r ) on the boundary. If so one could 
then also claim, in view of [23j . that the chiral conden- 
sate vanishes when averaged over longer distances. In the 
Skyrme model, approximate chiral symmetry restoration 
is relation to breakup of Skyrmions in half-Skyrmions, 
and it has been suggested in [55 that this breakup can be 
translated into a holographic mechanism for chiral sym- 
metry restoration. The precise implementation of this 
idea in our setup is however not clear to us. 



H. More than one flavor 

The Nf — 1 case is clearly much simpler than the 
Nf > 1 case where one has to deal with the non-abelian 
world-volume theory. The Nf — 1 case is also more sin- 
gular because there are only point-like instantons in an 
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abelian theory which we modeled using charged parti- 
cles. For Nf > 1 one can describe instantons in terms 
of smooth gauge-field configurations, though exact multi- 
instanton solutions in non-abelian Born-Infcld theory are 
not known. Isolated instantons in the Nf > 1 the- 
ory carry an "orientation" in flavor space, so for widely 
spaced instantons one can roughly think of them in terms 
of charged particles that carry extra spin-like quantum 
numbers. Interactions between instantons crucially de- 
pend on their relative orientation [32j [56] , so a mean-field 
description would probably be much more complicated 
and one would need to keep track of the extra quantum 
numbers. One could try to get a very crude description 
of the Nf > 1 case by taking the Nf — 1 case and simply 
adjusting the mass and charge of the baryons to account 
for extra contributions to the interactions coming from 
the orientation. As we saw in section |Vj such a simple 
modification of the theory can already cause important 
qualitative differences. 



paired with the radial electric field induces an interesting 
mixing pattern between the axial and vector mesons car- 
rying non-zero momentum and living in a dense baryonic 
medium. An interesting spin-off of our work is that the 
masses of mesons might be significantly different from 
their vacuum values due to density enhanced interac- 
tions with baryons, which will affect at least the numer- 
ical findings of [35] when applied to the Sakai-Sugimoto 
model. 

Recently, an interesting attempt was made to evaluate 
the chiral condensate at non-zero baryon density [59] . 
We want to stress however that this paper focused on 
three-dimensional densities and it would interesting to 
generalize its findings to four-dimensional densities that 
we consider, as well as try to understand possible impli- 
cations of the granular nature of the holographic medium 
on the chiral condensate in the dual field theory at non- 
zero baryon density. 



I. Relation to other work on the subject 

Earlier works on dense cold Sakai-Sugimoto model in- 
clude [23 H3 1311 122] • The crucial difference between our 
approach and that of (26] [27] is that we take into ac- 
count the strong electromagnetic repulsion between the 
holographic baryons arising at C(A 2 ) densities and study 
four-dimensional baryon lattices along the lines of [3T] . It 
is fair to say that so far there is no satisfactory macro- 
scopic description of such lattices that starts from first 
principles. The only microscopic attempt that we are 
aware of is that of |32| . 

The Chern-Simons term-induced instabilities in holo- 
graphic system are a subject of significant recent research 
interests, but in the case of holographic QCD there are 
only a few such papers 28 - 30 , 33J EH] • The only article 
among these that considers the baryonic phase at zero 
temperature in the Sakai-Sugimoto model is [49], but it 
focuses on three-dimensional baryon lattices and does not 
take into account the necessity of checking normalizabil- 
ity of gauge field perturbations at both ends of the flavor 
brane that we discussed in section [V] 

To the best of our knowledge, the meson spectra at 
parametrically large densities were studied only in [57] 
(see also [57J [SS]). This paper also focused on three- 
dimensional densities of the holographic baryons and did 
not take into the direct holographic meson-baryon inter- 
actions into account that we discussed at length in section 
VI The effects taken into account in [27] were the nonlin- 
earities of the DBI action, but as shown in Fig. [l6]this is 
not the key mechanism that affects the masses of mesons. 
One needs to take into account though that the radial 
electric field generated in the three-dimensional approx- 
imation to the holographic density of baryons might be 
bigger than the one in our case for parametrically large 
densities. 

In 1331 it was noticed that the Chern-Simons term 



J. Concluding remarks 

We think that the most important outcome of our pa- 
per is the realization that parametrically large baryon 
densities, of order of the square of the 't Hooft coupling 
constant, may lead to interesting physical effects in the 
Sakai-Sugimoto model resembling to some extent the ef- 
fects expected to arise in the Quarkyonic phase. While 
some of these effects were already anticipated in the liter- 
ature, e.g. four-dimensional densities of the holographic 
baryons appeared already in |31j and were recently dis- 
cussed in |32j . some others are new, like the density en- 
hanced holographic meson-baryon interactions. The cru- 
cial obstacle that we encountered in our studies was the 
lack of controllable microscopic insight in the nature of 
the correct coarse-grained description. 

Although the main motivation of our work was to make 
contact with the conjectured Quarkyonic phase in large- 
N c QCD, we found neither the exact chiral symmetry 
restoration in the homogeneous phase nor the unstable 
modes breaking translational invariance. It might well 
be that these features arise once one uses the correct 
coarse-grained description. An indication of this might 
be the fact that tweaking the parameters of our effec- 
tive lagrangian for perturbations does in fact lead to the 
appearance of a particular unstable (actually marginally 
stable) modes. 

We think that the most interesting question raised by 
our studies is the mechanism of the possible chiral sym- 
metry restoration due to non-zero density. Hence, the 
most obvious open direction would be to perform a more 
detailed microscopic analysis of the baryon lattice, to un- 
derstand the fate of the pion at finite densities within 
these models, as well as to evaluate the chiral consensate 
more accurately. 
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Appendix A: Equations for meson perturbations 

To get the equations for mesonic perturbations we first consider generic perturbations 5A a x a = Sa 3 (z,t)dx 3 and 
expand the DBI Lagrangian to second order in the perturbation 



S meson = J dzd i x\^ 1 '\l + |!)Viy_ det( £ o6) _|9/2 (1+ ^^(^ M ab M ch F bc F ha ^, (Al) 

where M ab = (g a b + Fab)" 1 and F a b is the background field strength which has the only non-zero components 
F z q = —Fq z = dzAo(z). Using the Euler-Lagrange equations, direct coupling of meson perturbations to baryon 
distribution and the ansatz Sd 3 (z,t) — <5a 3 (z)e~ 4Wt we get the desired coefficients 
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+ 5zA' (zy + y'(z)(z(5 - 2£(z* + ?))j?{z) + 3(z 2 + + ? - l)y"{z))j 
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The density of baryons on the D8 branes takes a simple algebraic form and is given by 

4/3 
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where the coefficients are 
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D 2 (z) = 9ft(£ 2 + 5 2 ) (| 3 +|z 2 -l) (9(| 2 + i 2 ) (e 3 +|5 2 -l) 2 y'(z) 2 +4| 5 / 3 z 2 (| 2 + 5 2 ) 
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The determinant of the metric is 



D 3 (z) = -\^m~z (f + z 2 | 3 + £z 2 - 1 
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The embedding function, obtained from the equations of motion, is given by 



(A9) 
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(A10) 



As was explained in the text it is explicitly seen that the embedding function depends on the baryon density only 
through the background gauge field (and matching condition on the boundary of the charge distribution). In the 
antipodal case £ = 1 and y'(z) = which makes the shape of D8 independent of the baryon density. 
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